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ABSTRACT 

This is the second of two papers reporting observations and analysis of the 
unusually bright (m^ = 14.4), luminous {Mb = —25.5), nearby {z = 0.192) 
narrow-line quasar PHL 1811. The first paper reported that PHL 1811 is in- 
trinsically X-ray weak, and presented a spectral energy distribution (SED). Here 
we present HST STIS optical and UV spectra, and ground-based optical spectra. 
The optical and UV line emission is very unusual. There is no evidence for forbid- 
den or semiforbidden lines. The near-UV spectrum is dominated by very strong 
Fe II and Fe III, and unusual low-ionization lines such as Na I D and Ca II H&K 
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are observed. High- ionization lines are very weak; C IV has an equivalent width 
of 6.6A, a factor of ~ 5 smaller than measured from quasar composite spectra. 
An unusual feature near 1200A can be deblended in terms of Lya, N V, Si II and 
C III* using the blueshifted C IV profile as a template. Photoionization modeling 
shows that the unusual line emission can be explained qualitatively by the unusu- 
ally soft SED. Principally, a low gas temperature results in inefficient emission 
of coUisionally-excited lines, including the semiforbidden lines generally used as 
density diagnostics. The emission resembles that of high-density gas; in both 
cases this is a consequence of inefficient cooling. PHL 1811 is very unusual, but 
we note that quasar surveys are generally biased against finding similar objects. 

Subject headings: quasars: emission lines — quasars: individual (PHL 1811) 



1. Introduction 



Strong, broad emission lines are an identifying feature of Active Galactic Nuclei (AGN) 
optical and UV spectra. The gas emitting these fines is illuminated and photoionized by the 
continuum spectrum emitted by the accretion disk and corona in the central engine. The 
emission line spectrum in AGN is broadly similar from object to object over many decades 
of luminosity. This makes sense because photoionization and thermal equilibrium tends to 
produce strong recombination lines, as well as strong resonance lines from easily excited 
abundant ions. 

Different elements and ions in the line-emitting gas have a range of ionization potentials, 
so naively one would think that the variations in the spectral energy distribution should in- 
fluence the emission-line equivalent widths and emission-line ratios. The influence is muted 
by the fact that many of the metal ions are excited predominately by collisions, so that the 
ionization state of the gas depends on the ionization parameter U = ^/nc, where $ is the 
photon flux, n is the gas density, and c is the speed of light. This me a.ns that the depen- 



dence on spectral energy distribution should be a secondary effect (e.g.. iKrolik fc Kallman 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associated with proposal #9181. 

^Based on observations obtained at Kitt Peak National Observatory, a division of the National Optical 
Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under cooperative agreement with the National Science Foundation. 
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19881 : ICasebeer. Leighly &: Baron II2006I ). Second, the 



ine emitting region i s most probably 



extended and comprises gas with a range of densities ([Baldwin et al. Ill995l ). Another com- 
plication is that the broad-line region may not see the continuum that we see becau se some 
components of the central engine may radiate more isotropically than others (e.g., iNetzer 



1987f ). 



Despite these complications, evidence that differences in the spectral energy distribu- 
tion have an effect on the broad-line region emission and kinematics has been found recently. 
It is now accepted that the Baldwin effect (the empirical anti correlation betw een the BLR 
emission-line equivalent widths and the continuum luminosity; [Baldwin 1119771 ) is a result of 
the softening of the spectral ener gy distribution in more luminous objects as a result of a 



larger black hole mass (see, e.g., ICasebeer. Leighly &: Baron 1 12006| . for a review); the ob- 



served anticorrelation between the Ba ldwin effect slope and the ionization potential for each 
line is the most compelling evidence ( Dietrich et al. II2002I ). The accretion rate, which also 
may aff ect the spectral energy distribution, may also influence the emission-line equivalent 
widths ( iBaskin fc Laor II2004I ). and the fact that two properties may influence the equiv- 
alent width perhaps contributes to the large scatter in the Baldwin effect correlation. In 
addition, the s pectral energy distribution can affect the line fluxes and equivalent widths in 
indirect ways. iLeighly I (120041 ) discussed the fact that an X-ray deficient spectral energy dis- 
tribution better explained the broad, blueshifted high- ionization line ratios in two extreme 
Narrow-line Seyfert 1 Galaxies (NLSls). Then, she showed that the narrow, intermediate- 
ionization lines were better explained if the continuum incident upon them had been first 
transmitted through the gas emitting the high-ionization lines. In addition, they found 
that the strong high-ionization line emission, in particular, the strong O VI line, in the 
l ow-luminosity NLSl RE 1034-1-39 is a consequence of its hard spectral energy distribution 
( ICasebeer. Leighly fc Baron 1120061 ). 

This is the second of two papers that report optical, UV and X-ray observations and 
analysis of the nearb y {z = 0.192), luminous {Mb = —25.5) narrow-line quasar PHL 1811 
( ILeighly et al. 1 120071 ). P HL 1811 was first cata loged as a blue object in the Palomar- 
Haro-Luyten plate survey (IHaro &: Luyten Ill962l ). It then rediscovered in the optical fol- 
lowu p of the VLA Faint Images of the Radio Sky at T wenty Centimeters (FIRST) sur- 
vey (Iwhite et al. 1119971 : Becker. White, fc Helfand Ill995h . It is extremely bright (B=14.4, 
R=14.1); it is the second brightest quasar at z > 0.1 after 3C 273. Being so bright, it is a 
very good background source for studies of the intergalactic and interstellar medium; fur- 
thermore, a FUSE observ a tion found it to have a rare Lyman limit system which has been 
studied by IJenkins et al. I (120031 . 120051 ). It was odd, however, that such a bright quasar was 
not detected in the ROSAT All Sky Survey (RASS). In comparison with other quasars of 
its luminosity, the expected RASS count rate is about 0.5 s~^; we placed an upper limit of 
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1.3 X 10 ^ counts s ^ (ILeighly et al. Il200ll ). A pointed B eppoS AX ohservation detected the 



object, but it was still anomalously weak. Too few photons were obtained in the BeppoSAX 



observ ation to unambiguously determine the cause of the X-ray weakness; iLeighly et al. 



(I2OOII ) speculated that either it is intrinsically X-ray weak, or it is a nearby broad-absorption 
line quasar and the X-ray emission is absorbed, or it is highly variable, and we observed it 
both times in a low state. 



In the companion paper (ILeighly et al. 1120071 . hereafter Paper I), we reported the results 



of five X-ray observations of PHL 1811 using Chandra, XMM-Newton, and Swift. We also 
reported the simultaneous UV photometry using the OM on XMM-Newton and the UVOT on 
Swift, as well as optical photometry obtained at the MDM observatory. These observations 
confirmed that PHL 1811 is an anomalously X-ray weak quasar. The X-ray spectrum is steep, 
with photon index between 2 and 2.6, and there is no evidence for absorption in excess of the 
Galactic H I column. The X-ray flux varied by a factor of ~ 5 among the 5 observations, but 
the UV photometry, when compared with the HST spectrum presented and analyzed in this 
paper, showed no detectable variability. The inferred aox, defined as the point-to-point slope 
between 2500A and 2 keV, is measu red to be —2.3 ± 0. 1. Typical quasars with PHL 1811's 



optical luminosity have aox of —1.6 (ISteffen et al. 1120061 ). Accounting for scatter, this means 



that PHL 1811 is between 13 and 450 times fainter in X-rays than other quasars with the 
same UV luminosity. Along with the RO SAT All Sky Survey upper limit and the BeppoSAX 
observation, we have now observed PHL 1811 in X-rays seven times and it is always found 
to be X-ray weak. While we can never disprove the hypothesis that we coincidentally always 
observe it in a transient low state, the odds of that being true are decreasing. We concluded 
in Paper I that PHL 1811 is the best example of an intrinsically X-ray weak quasar. We 
present an updated spectral energy distribution of PHL 1811 in Paper I and discuss possible 
reasons for its X-ray weakness. 

In this paper, we present HST STIS UV and optical spectra, and ground-based optical 
spectra obtained at the KPNO 2.1 meter telescope (§2). As shown in §3, the emission-line 
properties of PHL 1811 are very unusual. There are no forbidden or semi-forbidden lines, 
and the high-ionization lines are weak. The spectra are dominated by low-ionization lines, 
and Fe II and Fe III are strong. In §4 we investigate, using Cloudy models, the physics of 
gas illuminated by an X-ray weak spectral energy distribution (SED). It turns out that gas 
illuminated by a soft SED has much different properties than gas illuminated by a normal 
SED, and furthermore, many of the peculiar features of the optical and UV spectra of 
PHL 1811, including the lack of semiforbidden lines, can be explained by this single factor. 
We summarize the principal results of the paper in §5. We also include a brief Appendix 
in which we explore additional properties of the soft SED th at are not directly related to 



the observations. Some of the results have been published in ILeighly. Halpern &: Jenkins 
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( 120041 ) ■ We assume a flat Universe with Ho = 70kms ^ Mpc ^ and Qyac = 0.73, unless 



otherwise specified. 



2. Optical and UV Observations and Analysis 

The HST observations were reduced using the standard pipeline. To reduce the effects 
of any possible fixed-pattern noise, four exposures with each of the chosen MAMA gratings 
were made with the target stepped along the slit. The separate exposures for each detector 
were each cross correlated with a preliminary average spectrum to search for any systematic 
offsets; none were found. The average spectrum was then computed using the IRAF task 
scombine and resampled to a linear binning. The HST STIS CCD observation was split into 
two separate exposures in order to aid in the rejection of cosmic rays. 



We presented an optical spectrum taken at the MMT telescope in 1997 in lLeighly et al. 



(120011 ). We obtained several other spectra at KPNO that have better resolution and signal- 
to- noise ratios that we present here (Table 1). Two of the KPNO spectra were taken within 
1.5 months of the HST and Chandra observations. 

To check the wavelength calibration of the UV spectra, we fit Galactic absorption 
lines (FUV: Si II A1260, A1527; C II A1335; C IV A1548, 1551; NUV: Mg II A2796, 2804; 
Fe II A2344, 2374, 2383, 2567, 2600). We find that the FUV spectrum is consistent with no 
systematic shift. However, all five lines in the NUV were offset from their rest wavelengths 
by 1.0-1.7 A. To rectify this, the wavelength scale in the NUV was adjusted blueward by 
1.35 A (the mean of the shifts of the seven NUV lines). 

There were no convenient absorption lines available to check the wavelength calibration 
in the HST STIS CCD and ground-based optical spectra. Assuming no systematic shifts 
does not lead to any anomalous or unexplainable results (see below) so we assume that the 
wavelength calibration in these spectra is satisfactory. 

The i/ST spectra are all given in the vacuum wavelengths; we convert the optical spectra 
to vacuum wavelengths uniformly for consistency. 

We next estimate the redshift of the object. There is no prominent emission feature in 
the FUV spectrum that could confidently be c onsidered to be emitte d at the rest wavelength. 



There is, however, a weak absorption line that iJenkins et al. I ( 120031 ) speculated originates in 



Lya in the rest frame of the object. This absorption line gives an estimate of the redshift of 
z = 0.1903. Interestingly, however, in the medium resolution i/ST spectra, there can be seen 
two lines near the Lya restframe; one is consistent with a redshift of 0.1901, and another. 
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smaller one that does not appear in the low resolution data is consistent with a redshift of 
0.1920. Alternatively, these may originate in galaxies other than the quasar host that are in 
the same group. 

In the NUV spectrum, there is a narrow emission line clearly originating in Fe II UV191. 
The three multiplet components are closely spaced, lying within 3A of one another; the gfij 
weighted average wavelength is 1786. 2A. Fitting this feature in the NUV spectrum gives 
an estimated redshift of z = 0.1900. In the HST CCD spectrum, we may consider Mg II 
to be emitted at the rest wavelength; it gives an estimate of the redshift of 0.1920. The 
Balmer lines in the optical spectra also give estimates near this value (0.1919 estimated from 
Ha in the July spectrum; 0.1924 estimate from H/? in the October blue spectrum; 0.1923 
from H/3 in the October red spectrum; 0.1921 from Ha in the October red spectrum). Note 
there is no observed [O III] or any other narrow-line region emission line. Thus, we see 
a difference in the estimated redshift between the UV and optical spectra by 0.002. This 
corresponds to wavelength shifts of 2.9A at observed frame Lya, 4.3A at Fe II UV191, and 
6.7A at Mg II. The estimated wavelength calibration uncertainty, obtained from the STIS 



data handbook flKim Ouijano et al. II2003I ). is 0.5-1.0 pixel in the MAMA and 0.2-0.5 pixel 
in the CCD. This corresponds to 0.3-0.6A in the FUV spectrum, 0.8-1.6A in the NUV 
spectrum, and 0.5-1.4A in the CCD spectrum - much smaller than the shift required to 
explain the Az = 0.002 between the redshift estimates. Therefore, we assume that the 
Balmer and Mg II lines are most likely to reflect rest wavelength of the object, we take the 
redshift to be z=0.192, and assume that the Fe II UV191 is blueshifted. 

We next construct a merged spectrum using the HST and ground-based spectra. For 
the i/ST spectra, we start at the far UV end and successively adjust the normalization and 
construct average, resampled spectra in the overlap regions. The absolute photometry is 
estimated to be 4% and 5% for the MAMA and CCD spectroscopy, according to the STIS 
data handbook. We found good agreement in the overlap regions throughout the FUV, NUV 
and optical HST spectra. 

Turning to the ground-based spectra, we found that the October red and blue KPNO 
spectra were consistent in slope and normalization. However, the July red spectrum had 
better SNR than the October red spectrum, and we wanted to use it. It had a slight 
difference in slope compared with the July spectrum. We resample the overlap region, which 
was most of both spectra, onto the same wavelength binning, divided the two, fit a quadratic 
function to the result and correct the slope difference in the July spectrum. We then merge 
the October blue spectrum with the corrected July red spectrum. 



Next, we merge the HST and KPNO spectra. The i/S" T spectrum is brighter by about 
30%, most likely because of slit losses in the KPNO spectra, and there was also a small 



- 7- 



difference in slope. We resampled the overlapping regions onto the same wavelength scale, 
fitted the ratio to a 9-node spline and corrected the KPNO spectrum to match the slope of 
the HST spectrum. 

Finally, we deredde n the merged merged spectruni usiri g E{B — V) = 0.046, estimated 



from t he infrared cirrus (jSchlegel. Finkbeiner. fc Davis Ill998l ) and the lCardelli. Clayton fc Mathis 
(119891 ) reddening law, and apply the redshift correction {z = 0.192), and 



The merged spectrum is shown in Fig. 1. Except for the small absorption line near rest- 
frame Lya, the other prominent a bsorption lines are all from our Galaxy or from intervening 
galaxies; they are all identified in iJenkins et al. I (120031 ) . and the reader is referred to that 
paper for the identifications. Emission lines expected from active galaxies are marked. We 
also mark emission lines from a particular type of ion - those in which the lowest excited 
level has A^* 7^ compared with the ground state, but rather that transition to the ground 
state is semiforbidden. These ions include Si"*"^, Si"*"^, C"*"^, and Al"*"^. The emission lines that 
we mark are those that have the first excited level as the lower level. This type of line is 
analogous to auroral lines, except that the metastable transitions to ground are semiforbidden 
rather than forbidden. The motivation for looking at these lines is discussed in §3.3.2. 



3. Emission-line Properties 

In this section, we discuss the optical and UV emission-line properties of PHL 1811. We 
organize the discussion by moving shortward in the spectrum: we first discuss the optical 
emission lines, then the near UV, and finally the far UV lines. 

To model the spectra, we use a combination of the IRAF task Specf it, and programs 
that we developed ourselves written in IDL. Specf it requires 1-sigma error bars on the 
spectra, which are not available for the optical spectra. The signal-to-noise ratio in the 
continuum of the optical spectra is greater than 30, so fit results are not dominated by 
statistics. We estimate uniform 1-sigma errors as the standard deviation in line-free segments 
of the data, and interpolate between these estimates. 



3.1. Optical Emission-line Properties 

We use the merged spectrum to analyze the emission-line properties. Fig. 2 displays 
the region of the spectrum near H/3 and Ha. The strong emission lines from Fe II that 
are common in NLSls are clearly seen. In order to determine which lines besides Fe II are 
present in the spectrum, and to measure the properties of the Balmer lines, we perform 



- 8 - 



an Fe II template subtraction, as is commonly done (e.g., iBoroson fc Green Ill992l : iLeighly 
1999bl ) using a template obtaine d from the prototyp i cal N LSl I Zw 1. We use two different 
templates. The first one is the iBoroson fc Green I (119921 ) template with a few additional 



iron lines added in the blue by D irk Grupe. The second one was recently developed by 
Veron- Getty. Joly. fc Veron I (120041 ). Each template is broadened by convolution with a 
Gaussian with FWHM = 400kms~^, scaled appropriately, and subtracted from the spectra. 
Neither template matches the Fe II emission in the region of H/3 very well (Fig. 2). We note 
that this might be expected for the Veron-Cetty, Joly & Veron template, since those authors 
remove forbidden Fe II emission, leaving only permitted Fe II in their template. 



The IBoroson fc Green I (119921 ) template produces a better fit in the vicinity of H/3, so 
we use that spectrum for further analysis in that region. I Veron- Getty. Joly. &: Veron I (120041 ) 
identified Fe II usually hidden in the base of Ha, and therefore their template yields an Ha 
line with a narrower, more symmetric base. The Ha and H/? line profiles have extended 
wings and clearly can not be modeled using a Gaussian. We use a Lorentzian profile plus a 
local linear continuum to fit them. The results are given in Table 2. 

We clearly observe H7, and report the results of the Lorentzian fit to this line in Table 
2. However, because the Fe II subtraction is really not very good in the region of this line, 
we do not consider these results to be very reliable. 

The optical spectra contain very few other lines besides hydrogen Balmer lines and Fe II. 
He II A4686 is not apparent, although it could be present if it is weak and masked by the 
imperfect Fe II subtraction. There is no trace of the generally strong narrow-line region line, 
[O III] A5007, or any other narrow-line region line. There is an excess near 6350 A that can 
be seen after subtraction with the Veron-Getty, Joly, & Veron template (Fig. 2). Based on 
the line identifications for I Zw 1 made in Veron-Getty, Joly, & Veron (2004), this could be 
Si II AA6347.10, 6371.40. 

The most interesting lines that are seen in the optical part of the spectrum are the 
very low ionization lines Na I D and Ga II H& K The Na I D AA5 889.89, 5895.89 doublet is 
frequently blended with the He I A5875.70 (e.g., [Thompson Ill99ll ). but in this case, as seen 
in Fig. 3, we are fairly confident that the entire feature is Na I D and no He I is present. We 
modeled the feature in two ways, using the Veron-template-subtracted spectrum. First, we 
modelled the doublet with two Lorentzians of equal flux, width and fixed separation. The 
results of this model comprises the entry for Na ID in Table 2. Notably, we find that the 
mean wavelength of the doublet, which would be the weighted mean of the feature if the gas 
were very optically thick to the hne, is equal to the lab value. The second model includes 
an additional Lorentzian to model He I also. The separations of the lines were fixed, and 
the widths and fluxes of the Na ID lines were constrained to be equal (as before). Since we 
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do not see He I, we constrain the fit by assuming that the width is the same as that of Ha, 



altho ugh we note that it is observed to be somewhat broader in Seyfert galaxies (j Crenshaw 
19861 ). The improvement in the fit was not significant (Ax^ = 1.2 for 201 fitted points), and 



the uncertainty on the fiux of the He I component is almost as large as the fiux (i.e, relative 
uncertainty is 90%). In addition, the mean wavelength for the Na ID line is now redward 
of the lab wavelength, although it is still consistent with that value. We conclude that He I 
is not detected, and evaluate the upper limit by varying the normalization until the is 
larger by 6.63 (99% confidence for one parameter of interest). In addition, we search for 
He I A7067. We find no evidence for it, and we estimate the upper limit in the same way. 
The results are given in Table 2. 



The weakness of He I A5876 is physically interesting. I Crenshaw I (119861 ) t abulate the 



fluxes of He I and Ha for 9 Seyfert galaxies; the mean He I/Ha ratio is 0.061. [Thompson 



( 1l99ll ) tabulate He I/Ha ratio for 6 quasars, and the mean ratio is 0.025. We measure an 
upper limit on the ratio of 0.0062 in PHL 1811. Such a low ratio is difficult to explain using 
photoionization models; we address this issue in §4.4. 



The strength of Na ID is also interesting. [Thompson I (jl99ll ) discussed the presence 
of this line in quasars. The ionization potential of Na"*" is only 5.14 eV, so rather special 
conditions are necessary to produce this line. Using Cloudy models, he found that a very 
high column density was necessary to shield the neutral sodium from photoionization in 
order to produce sufficiently large Na ID/Ha. 

Turning to the Ca II H& K doublet, we initially modelled these lines as Lorentzians, as 
the other lines are. The line widths are then found to be only 770 ± 83kms~^, about half 
the width of the other lines. Gaussian profiles fit these lines somewhat better (Ax^ = 15.4 
for 152 points), and the widths are still somewhat narrow (980 ± 80kms~^). Indeed, if we 
overlay their profiles on that of B.(3 (Fig. 4), they appear to be narrower than that line. The 
Na I D lines are also somewhat narrower than Ha (1365 ± 180kms~^), although we model 
those with a Lorentzian profile. 



Observations of the Ca II infrared triplet in AGN were reported by iPersson I (119881 ). 
He found that these lines were detected in at least nine of the forty objects surveyed, and 



that the line strengths were roughly correlated wi th the optical Fe 



ports of observations of Ca II H&K are rather rare (IFerland fc Persson 



I. Int erestingly, re- 
Il989h. although we 



note the recent detection of the se lines in IR AS 07598+6508 by IVeron-Cetty et al. I (120061 ). 
Ferland fc Persson I (119891 ) and jjoly I (11989! ) present models that show that very large col- 
umn densities are required to produce the observed infrared triplet emission. There are two 
hypotheses for the weakness of the Ca II H&K lines in most AGNs. First, the structure 
of the Ca"*" ion is similar to the Fe"*" ion. In each case, there are resonance transitions to a 
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high level, which decay back to the ground state, or to an intermediate metastable state. 
Some of the UV Fe II lines and the Ca II H&K lines correspond to the resonance transitions 
to the high level, while the optical Fe II lines and the Ca II IR triplet correspond to the 
decay to the intermediate metastable state. Thus, just as scattering of th e UV Fe II lines 
reprocess the iron emission into optical Fe II lines (IKwan fc Krolik Ill98ll ). scattering in a 
high-column-density gas converts the Ca II H&K lines into the infrared triplet. Alternatively 
or in addition, absorption in the host galaxy may make Ca II H&K in emission difficult to 
detect. Regardless, the similarities between Ca"*" and Fe"*" suggest that the study of Ca II may 
be a profitable way to understand the physical conditions of the low-ionization line-emitting 
region in AGN, as Ca"*" is a much simpler atom than Fe"^. 

The Ca II lines are well resolved, so in principle we can try to constrain the optical 
depth of the line-emitting region. We fit the lines with equal width and a fixed separation, 
but allow their fiuxes to be independent. We find that the K component fiux is then larger 
than the H component fiux, and the ratio between the two is measured to be 1.28. The gf 
values are 1.36 and 0.66 for K and H, respectively, thus predicting a ratio equal to 2 for 
optically thin gas. We note that there is some uncertainty in the measurement of the line 
fiuxes because the continuum is not smooth and is difficult to identify; however, the line 
ratio can be seen by eye to closer to one than to two. iLaor et al. I (Il997bl ) measure the ratio 
between the similar Mg II doublet components in I Zw 1 to be 1.2. Mg II is isoelectronic 
with Ca II and these lines are both transitions from the lowest-level excited state to the 
ground state, so the interpretation of the low ratio should be similar. In both cases it means 
that that the gas is sufficiently optically thick that enough scatterings occur to redistribute 
the lines; the ratio would be 1 to 1 if completely optically thick. 

The ratio of Fe II to H/? is a parameter frequently u sed t o describ e d the s trength of the 
Fe II emission. As was done in iBoroson & Green I (119921 ) and iLeighly I (jl999bl ). we integrate 
the subtracted Fe II model between 4344 and 4684 A. Dividing the fiux measured from 
the model by the integrated and modeled H/3 fiuxes yields an es timated ratio in the range 
of 1.22-1.35. Comparing with the sample of NLSls reported in iLeighly I (jl999bl ). we find 
that the estimated ratio for PHL 1811 is in the middle of the range reported there, but 
somewhat larger than the median and mean (1.0 and 1.1, respectively). Thus, PHL 1811 is 
a moderately strong Fe II emitter. 

In summary, we find that the optical spectrum from PHL 1811 contains prominent 
Balmer and Fe II lines, similar to many other NLSls. The spectrum is somewhat unusual 
in that there are no forbidden lines; in particular the generally strong [O III] A5007 line 
is absent. NLSls are identified by their weak forbidden lines, and a few other objects 
appear to have no [O III] also (e.g., RX J0134.2-4258, ICrupe et al.1l200oh . There are 



no high-ionization coronal lines which sometimes appear in optical spectra of NLSls (e.g., 
Dietrich. Crenshaw. &: Kraemer 1 12005| ). PHL 1811 is unusual in that it shows prominent 
very low-ionization lines from Na ID and Ca II H&K, but no He I. 



3.2. The Near-UV Spectrum 

The near UV portion of the spectrum of AGN, between around 1600 A and 3500 A, 
generally exhibits strong and prominent emission lines from C III] A 1909 and Mg II A2800. A 
pseudocontinuum of Fe II may also be present over the entire region, and is most prominent 
between 2200 and 3050A. Fe III and Fe I may also contribute. Good signal-to-noise spectra 
often reveal contributions from numerous weak semiforbidden lines of oxygen, nitrogen, 
silicon, and carbon. The near UV spectrum from PHL 1811 is different from many AGN in 
that no semiforbidden or forbidden lines are observed. The spectrum is strongly dominated 
by Fe II and Fe III. 



3.2.1. Mg II, Fe II and Fe III in PHL 1811 



The properties of UV Fe II and Mg II in Active Galactic Nuclei (AGN) are important 
diagnostics f or several reasons. Lumino us quasars are probes of the early Universe. As 
discussed by iHamann fc Ferland I (119931 ) and others, the production of iron is thought to 
lag that of the a elements including magnesium due to different formation mechanisms: 
magnesium and about half of the iron (INomoto. Nakamura. Kobayashi Ill992l ) are produced 
in supernovae of massive, rapidly-evolving stars, while the remainder of the iron is produced 
largely in Type la supernova which are triggered by accretion onto a white dwarf, and are 
delayed with respect to core-collapse supernova. Observation of an evolution of the [Fe/Mg] 
ratio with redshift could constrain the time of the first burst of star formation in the Universe. 
The UV Fe II and Mg II are found conveniently near one another in the rest-UV bandpass, 
and they should both be produced in gas with similar properties, so it was thought that the 
Fe II/Mg II ratio could be an abunda nce diagnostic. No cle ar evidence for Fe II/Mg II ratio 
evolution has been observed yet (e.g.. iDietrich et al. 1120031 ). 



As in the optical ba ndpass, the flux and equiyalent width of Fe II are frequently measured - 



using a template (e.g.. ICorbin &: Boroson I Il996l : iForster et al. 1 1200 ll : iDietrich et al. 1 12002 



Leighly &: Moore 1 120061 ). We developed a UV Fe II te mplate from the HST spectru m of 
the prototypical Narrow-line Quasar 1 Zw 1, following IVestergaard &: Wilkes I (l200ll ): we 
subtracted a power law identified at the relatively line-free regions near 2200A and 3050A: 
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absorption lines and prominent emission lines not attributable to Fe II were then subtracted. 
Using this template, we recently measured the Fe II/Mg II ratio in the spectra from 903 



interm ediate-redshift narrow-line quasars from the Sloan Digital Sky Survey (ILeighly fc Moore 



20061 ). In most of those spectra, the template modelled the iron quite well in the 2200-3090 
A bandpass. However, when we try to apply the template to PHL 1811, we cannot obtain 
a good fit. While the region of the spectrum around Mg II matches the template well, there 
is significant excess emission at shorter wavelengths. To illustrate this problem, we adopt 
a strategy frequently used by X-ray astronomers: we fit the spectrum with a linear model, 
the iron template, and two Gaussians to model the Mg II doublet in the ~ 2650-3090 A 
bandpass. Then, we extrapolate the model to shorter wavelengths. The result is shown 
in Fig. 5a for PHL 1811. It is interesting that the template and continuum match the dip 
around 2250A fairly well, but there is significant excess emission between 2300A and 2500A 
and shortward of 2250A. 

For comparison, we perform the same analysis on the spectra of several other quasars and 
a composite spectrum. In Fig. 5b we show I Zw 1. Naturally, the t emplate fits the spectrum 



perfectly. Next, we show the largest Fe II/Mg II composite from iLeighly fc Moore I (12006 
shown in their Fig. 6). In this case and for the next case, we ignore the wavelengths containing 
the prominent emission from Fe II UV 62 and UV 63, low-lying transitions that sometimes 
produce strong emission near 2760 A, and trim the fitted bandpass. Finally, we show the 
spectrum of SDSS J105023. 68-010555. 5. The template matches the data rather well for both 
of these spectra. In Fig. 5a, however, we compa re PHL 1811 with t he HST FOS spectrum 



of the luminous NLSl RX J0134.2-4258 (e.g., iGrupe et al. Il2000f ). and the spectrum of 
SDSS J094257.80-004705.2. In both of these objects, like PHL 1811, extrapolation of the 
template leaves prominent excess emission at short wavelengths. 

What is the origin of the short wavelength excess in PHL 1811 and the other objects? In 
Fig. 6, we show the spectrum of PHL 1811 along with the upper level energies of transitions 
of Fe II and Fe III between 1000 and 3500 A; the atomic data was obtained from Kurucz 
Atomic Line Databasfl Around 2000A, there are many Fe III transitions with large Enrstem 
Aij coefficients. These include the multiplets from low-lying levels such as Fe III UV 34 near 
1900A and Fe III UV 48 near 2060A, which are seen as excesses in the spectrum. But there 
are also many more transitions with upper level energies between 15 and 19 eV. Thus, it 
may be that strong Fe III contributes to the excess shortward of 2250A. 

Between 2300-2500A, there are no high Aij transitions of Fe III. Fe II is more prominent 
in this region. The principal low-lying transitions with high A^j values between 2300 and 2500 



3 



http: / / cfa- www. harvard. edu/ amdata/ampdata/kurucz23 / sekur.html 
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A are UV 2 near 2400 A, UV 3 near 2340 A, and UV 35 and UV 36 near 2380 A. With upper 
level energy of only 5-6 eV, these are candidates for the short wavelength excess emission. 
However, UV 1 near 2610 A, UV 62 and 63 near 2750 A, and UV 64 near 2570 A also have 
upper level energies near 5-6 eV, and so one might expect these to be enhanced also. These 
transitions, and especially UV 62 and UV 63 are sometim es quite strong in some narrow-line 
quasars (e.g., SDSS J105023.68-010555.5 in Fig. 5b; also EaLeLaDEoi), but they do not 
appear to be excessively strong in PHL 1811 or in the objects in Fig. 5a. However, Fig. 6 
shows that there exists many higher-excitation Fe II transitions with high Aij values that 
produce emission in this region. Especially near 14 eV, there is a large concentration that is 
centered between 2300-2500 A, and we propose that it is this high-excitation Fe II emission 
that contributes to the Fe II pseudocontinuum so strongly in PHL 1811 and similar objects. 

Combined with the strong Fe III emission, the strong Fe II emission from high energy 
levels suggests that the iron in PHL 1811 and similar objects is characterized by a relatively 
higher ionization/excitation than in other objects. At first sight, this is a puzzling result, 
because the soft spectrum implies that the energy per ionizing photon is relatively low; using 
the PHL 1811 continuum developed for Cloudy simulations, we obtain an average energy per 
photon of only 19 eV. Also, the lack of hard X-rays implies that the partially ionized zone, 
where Fe II emission is thought generally to be emitted, will be quite shallow. But as we 
show in §4.3, it turns out that strong Fe II and Fe III are consequence of the soft spectral 
energy distribution through several mechanisms: the gas in the H II region is dominated by 
intermediate-ionization species, including Fe"*"^, there is strong continuum pumping of the 
Fe"*" ion, and there is strong Lya pumping associated with strong continuum emission in the 
partially-ionized zone. 

The Mg II is rather weak in PHL 1811, although in some sense it appears strong because 
it is prominent in comparison with the other lines. We find that the equivale nt width is just 



12.9 A. In c omparison, the equi valent widt hs from composite spectr a are 50A (IFrancis et al. 



19911), 64A f lZheng et ah Ill997h . and 34A flBrotherton et al. 1120011 ). In §4.1.3 we show that 
the low equivalent width for this collisionally excited line is plausibly a consequence of the 
low temperature in the line-emitting gas. 

Finally, we measure the UV Fe II/Mg II flux ratio, which, as noted before, is used as 
an abundance diagnostic in quasars. As shown above, the I Zw 1 template does not fit the 
UV Fe II spectrum in PHL 1811 well. Therefore, to measure the flux of this component, we 
define the continuum, subtract a pow erlaw model, integra te over the remainder from 2000 to 
3000A (the Fe II UV band defined by IVerner et al~l |20 0^, and subtra ct the flux of the Mg II 
line. Defining the continuum is difficult. As noted by IVerner et al. I (120041 ). there is a gap 
in the Fe II emission at around 3050A, but then a pseudocontinuum of varying amplitude 
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continues to at least ~ lOOOA. That pseudocontinuum is relatively low at 1800 or 2000 A, 
so to obtain a lower limit on the Fe II emission, we define two continua: from 1800-3050A 
and from 2000-3050A. The Fe II UV/Mg II fiux ratios obtained this way are 14.9 and 13.7 
for the 1800-3050 and 2000-3050A continua, respectively. These estimates are uncertain 
because some Fe III emission may be contained in the Fe II fiux estimate, which would tend 
to decrease the ratio. In addition, there is some uncertainty in the Mg II emission-line fiux, 
since Fe II template that was used assumes no Fe II emission un der the Mg 11; if there w ere 
Fe II under Mg II, the Fe II/Mg II fiux ratio would be larger (see lLeighly &: Moore II2006I . for 
a discussion of this point). For comparison, we found a maximum likelihood mean and la 



di spersion of the Fe I 



m 



/Mg II ratio of 3.93 ± 0.95 in the 903 SDSS quasar spectra analyzed 



Leighly &: Moore I (120061 ). Note that in that paper the Fe II fiux was integrated between 



2200 and 3050A. Thus, PHL 1811 shows an enhanced Fe II/Mg II fiux ratio, and enhanced 
UV Fe II emission compared with other quasars. 



3.2.2. No Semi- forbidden Lines in PHL 1811 



Fig. 7 shows a portion of the near-UV spectrum from PHL 1811. The HST spectrum 
of the prototypical NLSl I Zw 1 is also shown. The positions of permitted, semiforbidden, 
and forbidden emission lines commonly seen in AGN are marked. Prominent Fe II and Fe III 
lines are marked, although as discussed above, blends of lines from the same ions produce a 
strong pseudocontinuum. Absorption lines that are present all originate in our Galaxy. 

The most prominent difference between these two spectra is the relative weakness of 
all emission features in the spectrum of PHL 1811. In most AGN, C III] A1909 is one of 
the strongest emission lines; the equivalent width of the 1900 A feature in quasar composite 
spectra is around 20A (Francis et al. 1991; Brotherton et al. 2001). In I Zw 1, C III] is 
weak relative to S i III], perhaps because of a higher density in the broad-line region (e.g., 
Laor et al. I Il997bl ). In the spectrum of PHL 1811, there is a bump near the expected 
position of C III] . However, as the inset figure shows, the central wavelength of that feature 
is somewhat longward (~ 1915 A) of the expected wavelength for C III]. Its wavelength 
matches that of the Fe III A1914.1 emission hu e, part of the Fe III UV3 4 triplet (1894.8A, 
1914.1A, 1926.3A) more closely. As discussed in iLeighly fc Moore I (120041 ) , the optically thin 
ratios of the li nes in this triplet should be 1:1.4:1.8. However, in some objects in the sample 
considered by iLeighly fc Moore I ( 20041 ). the 1914A component was strongly enhanced. The 
upper level energy of this transition corresponds to 1214. 6A, which is just l.lA from Lya. 
Therefore, it is quite plausible that if microturb ulence is present in th e line-emitting gas. 



this line may be strongly pumped by Lya (also [Johansson et al. Il2000l ). We suspect that 
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this phenomenon is occurring in PHL 1811, so that the feature near 1914A is Fe III, and 
C III] is not present. 

It is not straightforward to obtain an upper hmit on C III] because of the comphcated 
iron emission in the region. We assign the conservative upper hmit recorded in Table 2 as 
follows: we assume that our identification of the 1915A feature is mistaken, and it is really 
C III]. We identify as pseudocontinuum on each side of that feature and integrate over it to 
obtain the flux. 

The lack of semiforbidden lines in PHL 1811 naturally suggests that the emitting gas 
has a high density. Of the six semiforbidden lines (O III] A1664, N III] A1750, Si III] A1892, 
C III] A1909, N II] A2142, and C II] A2326) in the region of the spectrum show n in Fig. 7, 



Si III A1892 has the highest critical density (3 x lO^^cm"^; iHamann et al. II2002I ). However, 
as we show in §4, density indicators break down for gas illuminated by a soft spectral energy 
distribution (SED). In fact, the lack of semiforbidden lines can be completely attributed to 
the soft SED. Semiforbidden lines are collisionally excited, and collisionally-excited lines are 
weak in gas illuminated by a soft SED as a consequence of a low gas temperature. 



3.2.3. Other N ear- UV Line Emission 



We clearly see the Fe II UV 191 feature in the spectrum. This triplet (1785.277, 1786.758, 
1788. 004A) originates in a high-excitation state; the upper level energy is 9.8 eV. This 
line is prominent in stellar spectra; the excitation mechanisr n is not known. Dielectroni c 
recombination has been suggested as a possible mechanism (j Johansson fc Hansen Ill988l ). 
That mechanism might be consistent with our hypothesis, presented in §3.2.1, that the near- 
UV line emission is dominated by Fe"*"^ and high-excitation Fe II. We model this line with 
three Gaussian compon ents fixed at their expected optically thin relative fluxes of 2.0:1.5:1.0 
( Johansson et al. Ill995l ). The best fitting gf weighted wavelength of the feature is 1783. 2A. 
Because this feature it so narrow, it can be seen to be clearly blueshifted (by 570 km s~^) from 
the expected gf weighted wavelength of 1786. 4A. A blueshift in a narrow, low-excitation line 
is not expected. Perhaps this is telling us something about the emission mechanism? 

We also identify Al III AA1854.7, 1862.8. We model this feature with two Gaussians 
with equal widths, fluxes and a fixed distance between. We measure a relatively large width 
(Table 2; 3050 km s~^), and large redshift (1340 km s~^). While the large widths and shifts 
may be real, it is also quite possible that the line is blended with other emission lines. An 
attractive possibility are t he "unexpected U V" lines at 1870A and 1873A resulting from 
Lya pumping of Fe II (e.g., iRudy et al. Il2000l ). We note that there seems to be a small blue 
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shift from the mark ed rest-wavelength position in I Zw 1 for both Fe II UV 191 and Al III 
Jlaor et al. Ill997bh . 

We also notice an excess in the spectrum at 1750A that could be the resonance transition 
of Ni II. In most AGN, that line would be hidden under the N III] multiplet; however, in this 
object, since other semiforbidden lines are conspicuously weak, there is no reason to suspect 
that N III] is present. 

As will be discussed in ^3.3. 2, we infer prominent Si II emission from multiplets near 
1263A and 1194A. As discussed in [Baldwin et al. I ( 19961 ). if these emission lines are produced 
by collisional excitation, we expect to see emission from lower-lying resonance transitions near 
1815A and 1531 A, neither of which are observed. The 1531 A line is difficult to identify, 
because it is hidden under C IV. We obtain an upper limit to the 18 15 A line by fitting with 
the C IV template, described in §3.3.1, fixed at the expected wavelength. We then vary the 
fiux of the template until we obtain larger by 6.63, which corresponds to 99% confidence 
for one parameter of interest. We note that our upper limit is conservative; because the 
spectrum is complicated in this region, we are forced to constrain the fit over only a narrow 
bandpass and therefore cannot constrain the continuum very well. 



3.3. The Far-UV Spectrum 



The far UV continuum and emission lines are both unusual in PHL 1811. Typical 
quasars show a break to a fiatter spectrum at ^ lOOOA, as is seen in a composite spectrum 



constructed from HST spectra of quasars (IZheng et al. 1119971 : iTelfer et al. 



20021 ). A break 



in this v icinity is also observ ed in the spectral energy distributions of individual low-redshift 
objects (IShang et al. 1120051 ). which means that the break cannot be a consequence of Lya 
forest absorption. Fig. 8 shows the broadban d continuum spectru m of PHL 1811 in com- 
parison with the Francis composite spectrum ( Francis et al. Ill99ll ). We see no evidence for 
a break down to 970A, the extent of the HST spectrum. The continuum is possibly a power 
law over the entire optical-UV bandpass, with a large contribution from Fe II and Fe III 
between 2000 and 3200A. Alternatively, a rollover might be present at ~ 2500A, a longer 
wavelength than usual. 

Lya is generally a strong and prominent line in AGN spectra. Since it is a recombination 
line of hydrogen, it is expected to be strong in photoionized gas. The equivalent width of 
Lya plus N V in composite spectra is foun d to be ~ 50 to ~ 90 A ( Francis et al. Ill991 



Zheng et al. 1 119971 : iBrotherton et al. I l200ll ). In addition, C IV is generally observed to 



be very strong as well. C IV is expected to be strong because C"^^ is easily excited by 
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coUisions. Composite spectra reveal a C IV" equivalent width of 33-59A ( Francis et al. Ill991 



Zheng et al. 1119971 : iBrotherton et al. Il200ll ). In PHL 1811, C IV and a feature near Lya are 



observed, but their equivalent widths are small (Table 2; see below). 

Although PHL 1811 has an unusual far-UV spectrum, it is not a unique object. Fig. 
9 shows the spect ra from two other high-luminosity narrow-line quasars: RX J0134.2— 4258 
(IGoodrich 1 120001 ) and PHL 1092 (Leighly et al. in prep.). All three objects show low 
equivalent - width emission lin es and blue far-UV spectra. For comparison, t he HST composite 
spect r um (IZheng et al. 1119971) . the spectrum from the NLSl IH 0707—495 (ILeighly fc Moore 



2004J : iLeigh 



1. 



200 



3), and a sp ectrum from the well-studied broad-line Seyfert 1 galaxy. 



NGC 5548 flKorista et al. Ill995f ) are shown. 



3.3.1. 



The CIV Line, He II A1640, and I4OO A Feature 



Fig. 10 shows the region of the spectrum that includes the C IV line. This line was 
somewhat difficult to analyze, because the continuum was difficult to identify. However, 
there is a relatively fiat region longward of the C IV line that stretches to 1800A that we 
could plausibly identify as the continuum. 



As discussed in lLeighly et al. I(l200ll ). and in Paper I, PHL 1811 is X-ray weak . iBrandt. Laor. fc Wills 



(I2OOOI ) investigated soft X-ray weak objects among the iBoroson fc Green I (119921 ) sample of 
PG quasars. They found that soft X-ray weak objects frequently show significant C IV ab- 
sorption lines; in fact, they showed that there is a trend for the most X-ray weak objects to 
have the largest equivalent- width absorption lines. They postulated that soft X-ray weakness 
is caused by photoelectric absorption of the X-rays by the same gas that is responsible for the 
C IV absorption lines. Indeed, in many cases, deep X-ray observations find t hat the signature 



of ab sorption is present in the X-ray spectra of soft X-ray weak objects (iGallagher et al. 



2OO2I ) . Since PHL 1811 is X-ray weak, we might expect to observe C IV absorption lines, 
and based on lBrandt. Laor. &: Wills I (l2000l ). we might expect an absorption line equivalent 
width of up to 20A. But as seen in Fig. 10, although the spectrum is noisy, no evidence 
for a C IV absorption line is seen. There are negative fiuctuations at the positions of reso- 
nance lines from gas in our Galaxy, namely Si II A1808 (seen at 1516. 8A) and possibly from 
Al III AA1854.7, 1862.8 (seen at 1556 and 1562.8 A). There is another negative fiuctuation at 
1529.3 A that we ide ntify as C IV absorption from the z{abs) = 0.17651 intervening system 
(j Jenkins et al. II2003I ). Thus PHL 1811 is not a typical soft X-ray weak object, and therefore 
we do not conclude that the soft X-ray weakness is a consequence of absorption. 



Fig. 10 shows the C IV emission line is blueshifted and asymmetric. Following the pro- 
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cedure discussed in iLeighly &: Moore I (120041 ). we create a template from the C IV line. 
The template fit, and the separate components, are shown overlaid on the C IV line. 
The shape of the PHL 1811 template is quite similar to the shape of the templates ob- 
tained IRAS 13224-3809 and IH 0707-495; given the lower equivalent width of the C IV 
line, and noiser spectra, they are essentially identical. To illustrate this, we overlay the 
IRAS 13224-3809 profile on Fig. 10. 

The equivalent width of the C IV l i ne is exceptionally low; we measure it to be just 
6.6A. As discussed in iLeighlv fc Moore I ( 20041). C IV equiva lent widths are much la r p;er in 
composite spectra; 



Francis et al. 



(ll99lh . lZheng et al. I (119971 ). and iBrotherton et al. I (l2001f ) 

measured C IV equivalent widths of 37, 59, and 33A respec tively. Fig. 8 comp ares the 
PHL 1811 C IV line with that from the Francis composite (jPrancis et al. 199ll ): it is a 
factor of > 5 sma l ler. T he two NLSls, IRAS 13224-3809 and IH 0707-495, analyzed in 
Leighly fc Moore I (120041 ) had low equivalent w idth of 13-15A; the C IV equivalent width in 
PHL 1811 is about half that. iLeighly I (120041 ) found that the weak C IV lines in those two 
NlSls, along with other unusual emission-line ratios, could be explained by a combination 
of high metallicity and an X-ray weak continuum. In §4 we show that the soft SED is 
responsible for the weak C IV in PHL 1811. 



P G 1407-1-265 is a quasar that has low-e quivalent width high-ionization lines (iMcDowell et al 



19951 ). iBlundell. Beasley. fc Bicknell I (120031 ) report the observation of a weak relativistically- 
beamed radio jet, and infer that the object is viewed face on. They then speculate that the 
low-equivalent width emission lines are a consequence of the face-on viewing angle, because 
they are diluted by the strong accretion disk emissi on. This reasoning h as been recently 
applied to another weak-line quasar, HE 0141—3932 (IReimers et al. II2005I ). 



While PHL 1811 may well be observed face-on, that cannot be the reason why the line 
equivalent widths are so low. First of all, one may no t be able to achi e ve a v ery large boost 
in continuum emission by viewing an object face on. iLaor fc Netzer I (119891 ) show that the 
ionizing flux from an accretion disk has nearly cos(6') dependence. Thus, assuming that the 
typical viewing angle to the accretion disk is ~ 30°, then one gains only a factor of 1.15 
by moving to a normal viewing angle. This is far from sufficient to explain the factor of 
~ 5 difference between the equivalent widths of C IV in PHL 1811 and the quasar composite 
spectrum. Of course, this analysis assumes that the continuum is not relativistically beamed; 
there is no reason to expect that it is, because although it is a radio source, it is a very weak 
one, and it is not radio-loud. 

Another argument against continuum boosting being responsible for the low high- 
ionization line equivalent widths is illustrated in Fig. 8. This figure shows that PHL 1811 
has a continuum slope quite similar to that of an average quasar; simultaneously the H/3 line 
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has a normal equivalent width and the C IV has a very low equivalent width. The reasoning 
goes as follows. First, if the continuum were uniformly boosted, one would expect both lines 
to have low equivalent widths; yet we see that H/3 has a normal equivalent width, while the 
C IV equivalent width is low. Second, if the continuum were boosted preferentially in the 
blue (e.g., a super-Eddington disk may allow r adiation to escape more easily along optically 



thin pores or channels (e.g., iBegelman 1 12002| ) in the hotter regions of the disk), then one 



would expect a bluer continuum than the average quasar. Instead, we see that PHL 1811 
has the continuum of a normal quasar, at least between C IV and H/?. 

We do not see any evidence for He II A1640. To estimate the upper limit of this line, 
we fit a C IV template shifted to the appropriate wavelength for He II. The best fit does 
not detect this line, so we estimate the upper limit by increasing the normalization until the 
has increased by 6.63 compared with the best fit (99% confidence for one parameter of 
interest). 

The 1400A feature is difficult to analyze. First, the continuum slope appears to change 
at ~ 1470A, although the continuum blueward of the 1400A feature appears to be flat and 
to extrapolate to the continuum shortward of the Lya feature. Second, the signal-to-noise 
ratio in that part of the spectrum is not very good. Finally, absorptio n lines from Galactic 



Al II A1670.79 and from C IV in the z=0.08093 damped Lya system f lJenkins et al. 1120031 ) 
cut through the middle of the feature. This feature is difficult to model, regardless, because 
it is usually a blend of the two multiplet lines of Si IV and the five multiplet lines of O IV] . 
Here, based on the lack of semiforbidden lines in the r est of the spectr um, it might be safe 



to assume that the entire feature is Si IV. However, in iLeighly I (120041 ) it was inferred that 
O IV] participated in the wind in IRAS 13224-3809 and IH 0707-495, and therefore was 
emitted in a different component than the intermediate-ionization lines such as C III] that 
are so clearly absent in PHL 1811. To emphasize the uncertainty in the identification of 
this feature, we refer to it as "unidentified" in Table 2. To estimate the fiux and equivalent 
width, we fit the entire feature with a Gaussian. The centroid is measured to be 1395 A, 
which possibly suggests an origin in Si IV, which has a gf weighted average wavelength of 
1396.75 A. 



3. 3. 2. The Lya /N V Region 



The far-UV spectrum of P HL 1811 is cluttered with absorption lines that originate in 
intervening absorption systems (j Jenkins et al. II2003I ). In order to better see the shape o f the 
far UV spectrum, we model out the absorption lines first. Using iJenkins et al. I (120031 ) as a 
guide, we fit local continua around each absorption line, and then model each absorption line 
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as a Gaussian. A comparison of the portion of the spectrum around Lya with and without 
the absorption hues removed is seen in the top panel of Fig. 11. 

The Lya/N V feature in PHL 1811 appears as a lump stretching from ~ 1140 to 
~ 1270A superimposed on a rising continuum that appears to be linear or a power law in 
the range 1090-1140 A and 1320-1360 A. The equivalent width of the whole feature is only 
24A, which is small compared with the 50-90A measurements of Lya+N V from composite 
spectra mentioned above. 

It is not obvious at first glance what the various bumps in the Lya/N V feature 
are. There are peaks that could correspond to Si II multiplets near 1193, 1260 and 1304 
A, and peaks tha t may be N V and Lya. We apply the approach previously used by 



Leighly &: Moore I (120041 ). and fit these features using the C IV template profile described in 
§3.3.1. A difference in the analysis presented here is that we now use the IRAF Specf it pro- 
gram which can fit an input profile, and properly minimizes to obtain the normalizations 
and spectral shifts. 

We first determine the continuum by fitting the line- free bands 1086-1103, 1116-1137 
and 1314-1350 A with a linear model. Initially, we assumed that the feature is composed 
of the five lines described above, and each one has the C IV template profile. We construct 
profiles for each multiplet, assuming first that the gas is optically thin and the multiplet ratios 
are proportional to the gf values, which were obtained from the NIST databas^. Initially, 
we fix the wavelengths of the profiles at the rest wavelengths, this implicitly assuming that 
gas emitting all the lines has uniform kinematics. The resulting fit is shown in the middle 
panel of Fig. 11. Many of the bumps are modelled, but not well. Specifically, the N V 
emission appears to be insufficiently blueshifted, and there is extra emission blueward of 
Si II A 1193. Emi ssion in this region has been previously identified as C III*A1175 (e.g.. 



Laor et al. Ill997bl ). so we assume that as the tentative identification of this excess. 



We also try profiles with all of the lines in the multiplets having equal strength, as would 
be appropriate if the gas were very optically thick. This model yielded a slightly larger 
but was not clearly distinguishable from the optically-thin case. We therefore proceeded to 
use the optically-thin profiles. 

We next apply the same model but add another template component around 1175 A 
to model the suspected C III* emission, and allow the central wavelengths of the profile 
components for each line to be free in the fitting. We obtain a much better fit, as seen in the 
lower panel of Fig. 11. The wavelengths, fiuxes and equivalent widths of this model are given 



"'http: // physics .nist .gov / cgi-bin/ AtData/ main_asd 
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in Table 2. We find an excellent correspondence between the laboratory wavelength and the 
observed wavelength of the template fits of Si II A1193, Lya and C III*; the remaining lines 
are somewhat blueshifted. The blueshift on Si II A1308 is conceivably partially due to an 
incomplete or faulty subtraction of a prominent absorption line. The blueshift on N V is 
plausibly real, as it might be expected that higher ionization lines may be emitted by gas 
with a higher velocity; this may be evidence that our template model is too simple. 

The presence of C III*A1175 is interesting. This line is e mitted as a transition to a 
metastable state that is the upper level for the C III] A1909 line. iLaor et al. I (Il997bl ) discuss 
possible reasons for a strong C III* line in the prototype Narr ow- line Seyfert 1 galaxy I Zw 1; 
we also tentatively identified this line in two other NLSls (ILeighly &: Moore 1120041 ). The 
presence of C III* suggests that there might be a similar transition from the metastable 
upper level of Si III] A1892 line, since Si"*"^ and C"^^ are isoelectronic. That transition should 
appear in the region 1108. 4-1113. 2A (the gf-weighted average wavelength is 1111. 59A). 
Indeed, there is a bump clearly present at these wavelengths that we confidently identify as 
Si III*A1112. Interestingly, it is not well fit by the template profile; we model it using first 
a single Gaussian, and then by the six components of the multiplet assuming optically thin 
ratios. For the second model, the width is 690 ± 130 km s (Table 2), and it is marginally- 
significantly blueshifted by 290kms~^. 

The successful deconvolution of the Lya/N V lump into these six components using 
the C IV profile as a template may be fortuitous, and we caution that it is certainly not 
unique and may not be complete. The presence of the metastable transitions of C'^^ and 
Si"*"^ led us to search for other similar metastable lines. On Fig. 1, we mark the positions of 
lines that have lower levels that are metastable, so that the transitions to ground state are 
semif or bidden, as described in §2. In the region of the Lya/N V lump, there are a couple of 
Si II* transitions that could be strong. We also mark Si III A1206.5, a resonance line that 
may also contribute; this line is analogous to C III A977. 



Strong high-excitation Si II has been previously reported by [Baldwin et al. I (119951 ). In 
§3.2.3, we obtained an upper limit on Si II A1814, and using the results of the deconvolution 
presented in this section, we measure a lower limi t on the Si II A1263/ Si II A1814 of 8.5, 
even larger than the v alue of 6 that was obtained by [Baldwin et al. I (Il995l ) from Q0207— 398. 
Baldwin et al. I (119951 ) showed that such high ratios are hard to explain, if the gas is collision- 
ally excited, as one expects that the lower excitation line to be stronger. They suggested that 
the high-excitation transitions are selectively excited. The key seems t o be that far UV lines 
are predominately excited by continuum pumping, generally speaking. iBottorff et al. I (120001 ) 
explored the effect of microturbulence on line ratios in AGN, and showed that since pumped 
lines are enhanced by turbulence, large Si II A1263/Si II A1814 ratios can be produced if 
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the turbulent velocity is high. Finally, ICasebeer. Leighly &: Baron I (120061 ) investigated Si II 
emission as a function of the shape of the ionizing continuum in their Appendix. They note 
that high-excitation Si II emission is stronger for softer continua, and that the Cloudy output 
information indicates that for soft continua, continuum pumping is important. 

In §3.3.1, we reported an upper limit on He II A 1640, and comparing with N V measured 
from the deconvolution pre sented here, we obtain a lower limit on the N V/He II ratio of 
11.6. iFerland et al. I (119961 ) reported N V/He II > 10 for "Component B" (the outflowing 
component) in Q0207— 398. N V and He II are both high-ionization lines from ions that 
require similar ionization potentials (He"*"^: 54.4 eV; N+^: 77.5 eV); therefo re, they should 



be em itted by gas with similar properties. A thorough analysis carried out by lFerland et al. 



( 119961 ) showed that the only way that the high N V/He II ratio can be produced is if the 
metallicity is enhanced by a factor of 5. A large N V/He II ratio in PHL 1811 may likewise 
imply a high metallicity also, although we again note the uncertainty in the deconvolution. 



3.3.3. Other FUVlmes 

There are a few relatively narrow, weak lines in the FUV that we try to identify. We 
were able to identify N II A 1085, which is a resonance line of N II. We model it with a single 
Gaussian, and also with the six Gaussians that comprise the multiplet. The latter fit yields 
a rather narrow line width, only 370kms~^. We also see another narrow line near 1034A, 
although this line is uncertain because considerable reconstruction had to be performed in 
the vicinity of this line. We were unable to confidently identify this line. It seems to be too 
centered and too narrow to be associated with O VI A1032, 1038. 

C II A 1335 is commonly seen in AGN, but we see no hint of it although it occurs in 
a relatively line-free region of the spectrum. To obtain an upper limit on the flux of this 
line, we add a C IV template component for this line to the FUV model, and fit over the 
entire region from 1086-1365A; we also include a short strip of continuum between 1425 and 
1462A. No significant emission is detected, so we estimate the upper limit on the flux of C II 
by increasing the normalization until the is larger than the best fitting value by 6.63 (99% 
confidence for one parameter of interest). 



4. Discussion 



The spectra and analysis presented in §3 reveal the unusual emission-line properties of 
the narrow-line quasar PHL 1811. We briefly summarize the principal results here. PHL 1811 
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has no forbidden or semiforbidden lines, and unusual very low-ionization lines are present 
(e.g., Ca II and Na ID). He I A5876 is very weak, and a very low He I/Ha ratio is inferred. The 
Fe II and Fe III pseudo continuum dominates the near-UV part of the spectrum. However, 
while the Fe II in the majority of AGN is fit well by the template developed from the 
high signal-to-noise HST observation of I Zw 1, the pseudocontinuum in PHL 1811 has a 
demonstratively different shape that we show is likely to be a consequence of the presence 
of high-excitation Fe II and strong Fe III. Some of the other low-ionization lines are strong, 
including the high-excitation Si II lines near 1194A and 1263A, but other low-ionization 
permitted line emission (e.g., Mg II) is weak or absent (e.g., C II A1335). While the hydrogen 
Balmer lines appear to have rather typical equivalent widths, the equivalent widths of the 
high-ionization lines are small, and in addition, the high-ionization lines have a blue-shifted 
profile. Deconvolution of a feature near Lya in terms of a template developed from the 
blueshifted C IV line reveals strong N V and high-excitation Si II, as well as C III* All 76. 

The goal of this section is to interpret these results. But first, we summarize previous 
work that has addressed some of the results. We note that previous work has been mentioned 
throughout §3, but we summar i ze the major points here also. The production of Na ID has 
been discussed by [Thompson I (1199 ll ). who inferred that since the ionization potential of 
Na"*" is only 5.14 eV, this line must be produced very deep in the partially-ioniz ed region 



where neutral sodium can survive . The production of Ca II has been discussed by iPersson 



f ll988h : iFerland fc Persson I fll989h andlJolv 



correlated with the optical Fe II ( Ferland &: Persson 



we report here, are rarely seen (IFerland &: Persson 



(Il989h 



Observationally, the infrared triplet is 



1989), but t he Ca II HfcK, the hues that 



19891). Both IFerland fc Persson I (119891 ) 



and lJoly I (119891 ) found that a high column density i s necessary to pr o duce the infrared triplet. 
Strong high-excitation Si II was first discussed by [Baldwin et al. I ( 119961 ). who pointed out 
the difficulty of producing large fluxes of the high-excitation lines such as S i II A1263, but 



at the same time producing very low fluxes of the resonance line Si II A1814. iBottorff et al. 



( I2OOOI ) showed that this problem may be solved if the emitting gas is turbulent, as the high 
excitation Si II lines are predominately pumped by the continuum and will be selectively 
excited by turbulence. Blueshifted high-ionization lines have been observed in other NLSls 
such as IRAS 13224-3809 and IH 0707-495 fe.g.. Ileighlv fc Moore l[20o3 ). The C IV lines 
in these two objects also have relatively low e quivalent widths, especially compared with 
the N V and the 1400A feature. iLeighly I (120041 ) showed that the gas emitting these lines is 
optically thin to the continuum, and that the line ratios can be explained by a somewhat 
X-ray weak spectral energy distribution and elev ated gas metal l icity. We do not detect 
He II A 1640 and thus infer a large N V/He II ratio. IFerland et al. I (119961 ) discussed the large 
N V/He II ratio in Q0207— 398, and found that the only way that this ratio can be explained 
is if the metallicity were enhanced compared with solar by Z ~ 5Zq. 
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What is the origin of the unusual emission-hne properties in PHL 1811? PHL 1811 
has one property that is much different than an ordinary quasar: it is intrinsically X-ray 
weak (Paper I), and therefore has an unusually soft spectral energy distribution. In the 
next several subsections, we examine the physics of gas illuminated by a very soft spectral 
energy distribution. We find that gas illuminated by a very soft SED has much different 
physical properties than gas illuminated by a typical AGN SED, and that most of the unusual 
observational features of PHL 1811 can be explained qualitatively by the soft SED alone. 
The low He II/Ha ratio cannot be explained by the soft SED, but we find that it can be 
explained by a "filtered" continuum; we discuss this idea in §4.4. Finally, in §4.5 we discuss 
whether there may be other objects like PHL 1811, and how we might identify them. 



4.1. Physical Conditions in Gas Illuminated by a Soft SED 

In this section, we investigate the physical conditions and line emission from gas illu- 
minated by an X-ray weak spectral energy distribution. Our goal is physical insight rather 
than a detailed match to the line emission, so we begin by exploring the physics of gas illumi- 
nated by a range of spectral energy distributions for a representative value of the ionization 
parameter (logf/ = —1.5) and density {logn = 11). These correspond to a photoionizing 
flux of log$ = 20. We use a hydrogen column density of loglNn) = 24.5, which corresponds 
to a thickness more than 100 times the depth of the hydrogen ionization front. We now 
define some terms. We refer to the region between the illuminated surface of the gas and 
the hydrogen ionization front as the H II region, and the region at greater depths than the 
hydrogen ionization front as the partially-ionized zone. In addition, we discuss the ionization 
potential to create emitting ions, and we refer to this as "ionization potential" , whereas more 
correctly, it is the ionization potential of the next lower charge state. 



4- 1.1. The Continua 
We use a range of continua to explore the influe nce of the spectral energy di s tribu 



tion. We start with the ones previously presented in in ICasebeer. Leighly fc Baron I (120061 ) 
(hereafter referred to as the CLE continua). They are described in detail in the Appendix 
of that paper; briefly, they constructed a number of spectral energy distributions parame- 
terized by the energy in eV of the UV bump high-energy cutoff. They assumed that the 
temperature of the cutoff Tcut-, which corresponds roughly to the inner edge of the accretion 
disk, is related to the luminosity by L oc T~^. Then, to determine the X-ray power law 
normalization, they assumed that the UV monochromatic luminosity at 2500A is related to 
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the X-ray monochromatic luminosity at 2 keV by ag^, , with ag ^, dependence on UV luminos- 



ity given by the regression developed by IWilkes et al. I (Il994j ). Finally, the proportionality 



constants were determined usin g the measured T„i.t from coordinated lUE and ROSAT oh- 



servations of the quasar 3C 273 (IWalter et al. Ill994j ). For more details, see the Appendix of 



Casebeer. Leighly fc Baron ( 20061 ). 



We also use the PHL 1811 SED developed from the HST and the brighter of the two 
Chandra observations, supplemented by a FUSE point from a non-simultaneous observation, 
and using the Cloudy AGN continuum at higher and lower energies where we have less 
information (Fig. 12). We do not know the shape of the extreme UV, so we simply join the 
FUSE point and the Chandra spectrum using a power law. We also use the Cloudy AGN 
continuum with parameters as used by Korista et al. (1997); hereafter referred to as "K97", 
and also shown in Fig. 12, this is taken to be a typical AGN continuum. 

Fig. 13 gives some information about the properties of these continua. The CLB continua 
are coded on the x-axis by their cutoff temperature, and the other two files are named 
"PHL1811", and "K97" for the PHL 1811 and Korista et al. continua, respectively. The top 
panel shows the aox for the continua. As noted in Paper I, PHL 1811 is anomalously X-ray 
weak and has a very steep aox = —2.3 ± 0.1. In comparison, the aox for the K97 spectrum 
is -1.4. 

In the next panel, we plot the continuum flux at 2500A assuming typical values of 
the input ionization parameter (logf/ = —1.5) and density (logn = 11). These values 
correspond of course to the same value of the ionizing photon flux, but the continuum flux in 
the optical or UV is much different. Specifically, an X-ray weak continuum must be brighter 
in the optical and UV in order to produce the same ionizing flux, since the ionizing flux is 
integrated from 13.6 eV toward shorter wavelengths. This will influence the equivalent widths 
of the lines produced; specifically, for objects that have soft continua like PHL 1811, the lines 
will appear to have lower equivalent width for the same global covering factor compared with 
an ordinary AGN continuum. The PHL 1811 continuum is 3.5 times brighter than the K97 
continuum at 2500A. This is certainly one factor contributing to the low equivalent widths 
of the line emission observed in PHL 1811. 

In the bottom panel of Fig. 13, we show the Compton temperature, which is a measure of 
the mean photon energy of the continuum. The gas would reach the Compton temperature 
if it were in equilibrium with the radiation field. The Compton temperature is very low 
for the PHL 1811 continuum, around 5 x lO^K; this temperature is similar to the typical 
temperature of a photoionized gas. 
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4.1.2. Temperature and Hydrogen Ionization 

We first examine the temperature structure of tlie ionized gas (Fig. 14). Tlie drop 
in temperature at around 1-2x10^^ cm occurs at tlie liydrogen ionization front; tlie kink at 
shallower depths occurs at the helium (He"^^ to He"*") ionization front. As might be expected, 
the very soft continuum from PHL 1811, and accompanying low Compton temperature, 
produce a rather cool temperature in the photoionized gas. The temperature in the partially- 
ionized zone for the PHL 1811 is exceptionally low; we will return to this point later. In 
contrast, the harder continua produce significantly higher temperatures, especially in the 
H II region. 

Next, we examine the fraction of ionized hydrogen as a function of depth (Fig. 15). The 
hydrogen ionization front is slightly deeper for the harder SEDs, perhaps marking additional 
ionization associated with the presence of soft X-rays, higher Compton temperature and 
greater heating. The fraction of ionized hydrogen in the partially ionized zone follows the 
temperature for the harder SEDs; that is, at a particular depth, the ionization fraction is 
higher for harder SEDs. That breaks down at the softest SEDS, where the ionization fraction 
rises again. For the PHL 1811 continuum, the ionization fraction is higher than for many of 
the harder SEDs. 

To understand the origin of the high ionization fraction in the partially-ionized zone 
of the gas illuminated by the PHL 1811 continua, we review the physical mechanisms that 
prod uce ionized hydrogen in the partia l ly-ionized zone . There are two principal mechanisms 
(e.g., ICoUin-Souffrin fc Dumont Ill989l : iFerland Ill999l ). The first mechanism involves X-ray 



photons that photoionize heavy elements through inner shell interactions creating secondary 
nonthermal electrons that ionize hydrogen. These electrons also heat the gas. The second 
mechanism involves photoionization of hydrogen in excited states (e.g., n = 2) by the direct 
continuum and the diffuse continuum. The hydrogen is in an excited state due to a circular 
process. Recombining hydrogen produces Lya; the optical depth to Lya in the partially- 
ionized zone is very high, so the Lya photons scatter many times before escaping the gas. 
That scattering may leave some hydrogen atoms in n = 2 long enough that they can be 
photoionized by the Balmer continuum of the incident quasar continuum. 

To investigate the importance of the latter process as a function of spectral energy 
distribution, we examine the rate of photoionization from n = 2 per hydrogen atom in 
= 2 by solving the radiation transport equation at any depth in the gas. We compute 
the photoionization rate from the incident continuum, attenuated by e""^", where t^, can 
be obtained using the extinction coefficient which can be output from Cloudy. We also 
compute the photoionization rate from the diffuse continuum, where the diffuse continuum 
Jj, is obtained by integrating the source function; this can also be output from Cloudy. The 



-27- 



resulting rates are shown in Fig. 16. For comparison, we show the coUisional ionization rate. 

Fig. 16 shows first that coUisional ionization is practically neghgible in the partially 
ionized zone. This is not surprising, since the temperature of the gas is a fraction of an eV. Of 
the photoionization rates, the contribution from the direct continuum is more important than 
that from the diffuse continuum, but the diffuse continuum contributes most to the harder 

contimia, from, e.g., ~ 15% for the K97 contimnim to ~ 40-60% in the very hardest continua. 
Considering the contribution from the direct continua alone, we see that the photoionization 
rate from n = 2 is highest for the softest continua. This is simply a consequence of the fact 
that for a particular ionizing flux, the flux density in the optical and UV will be higher for 
the soft continua (Fig. 13; §4.1.1). By itself, this predicts that the softest continua would 
produce the highest ionization fractions in the partially ionized zone. This is not what we 
see, so the difference must lie in the number of hydrogen in n = 2. 

The density of hydrogen in n = 2 can be output from Cloudy also, and we show it in 
left panel of Fig. 17. We see indeed that in the partially-ionized zone, this density is higher 
in objects with the hardest SEDs. However, this difference is insufficient to explain the 
dependence of the hydrogen ionization on the SED. For example, at a depth of 5 x 10^^ cm, 
the hydrogen fraction in gas illuminated by the PHL 1811 continuum is almost the same as 
in gas illuminated by the kT = 320 cV SED. At that depth, the total photoionization rate 
is a factor of five larger for PHL 1811, but the density of hydrogen in n = 2 is only a factor 
of 1.5 larger for the kT = 320 eV SED. Thus, the balance of hydrogen ionizations in gas 
illuminated by the hardest continua come from the first process mentioned above, the initial 
inner-shell X-ray photoionization and resulting photoelectrons. 

Next, we investigate the origin of the hydrogen in n = 2. We plot the ratio of the exci- 
tation temperature based on the populations in n = 1 and n = 2 to the electron temperature 
(Fig. 17). The ratio exceeds 1 in the partially-ionized zone, except for the gas illuminated by 
the hardest continua, which means that the temperature is not sufficiently high to explain 
the population in n — 2. Therefore, the radiation field must be responsible, and indeed, 
scattering of trapped Lya mentioned above is likely to be the origin of hydrogen in n = 2 in 
gas illuminated by the PHL 1811 continuum. 

The surprising aspect of Fig. 17 is that the gas illuminated by the PHL 1811 continuum 
exhibits a very large ratio, indicating an exceptionally high ratio of hydrogen in n = 2 
compared with that expected from a thermal gas. This trend continues to even higher n; in 
Fig. 18 we show the ratio of the density of hydrogen in n = 5 from the Cloudy output to the 
ratio expected based on the density in n = 1 and the electron temperature; this parameter is 
proportional to the departure coefficient for n = 5. This shows that departure coefficient is 
higher in gas illuminated by the PHL 1811 continuum by nearly three orders of magnitude 
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compared with gas illuminated by the hardest continua. 

The origin of the excited-state hydrogen in PHL 1811 is the strong diffuse continuum. 
In Fig. 19 we show the integrated diffuse continuum (i.e., energy in the diffuse continuum) 
multiplied by the depth, as a function of depth. The energy in the diffuse continuum is much 
larger for gas illuminated by the PHL 1811 continuum. Note that this does not necessarily 
mean that wc expect to see a strong continuum emission from this object, as the diffuse 
continuum used in this calculation is the emission at that depth, and the radiative transfer 
to the surface has not been folded in. This fact also explains the low electron temperature 
in the partially-ionized zone; the energy is locked up in the radiation field and excited states 
of hydrogen rather than electron kinetic energy. Thus, in PHL 1811, this strong continuum 
scatters in the partially-ionized zone, causing the hydrogen to be in high-excitation states. 
We will see why the continuum is so strong in the next section. 

To summarize this section, we find that the hydrogen in the partially-ionized zone in gas 
illuminated by the PHL 1811 continuum is rather highly ionized. This is a consequence of 
the bright Balmer continuum for a given photoionizing flux for soft SEDs, and the relatively 
large density of hydrogen in n = 2, which is a consequence of strong diffuse continuum. As 
will be shown in the next section, the low temperature means that the gas cannot cool by 
emission of collisionally-excited lines, so the energy is deposited in strong continuum and 
excited states of hydrogen, resulting in the relatively low gas temperature. 

4.1.3. Metal Ions: General Considerations 

The conditions in the line emitting gas can be understood further by examining the 
behavior of the metal ions. In Fig. 20 we show the ratios of the column densities of the 
ions of interest for the PHL 1811 continuum to the K97 continuum as a function of the 
ionization potential to create each ion. These ratios show that for ionization potentials less 
than ~ 10 eV, the column densities are approximately the same; the high point is Ca+, in 
which the column is a factor of 2.3 times greater in PHL 1811. For ionization potentials 
greater than 20 eV, the column densities are considerably lower in the gas illuminated by 
the PHL 1811 continuum. This is not surprising because the PHL 1811 continuum lacks 
high energy photons required to create these ions. For ionization potentials between 10 and 
20 eV, the columns in PHL 1811 are higher than in K97. This also makes sense because 
these are the ions that occupy the H II region, instead of the high-ionization species that are 
usually there. 

One might expect that if the column density of an ion is larger, the line emission from 
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that ion must also be larger. Fig. 20 shows the ratios of predicted emission-line strengths 
from gas illuminated by the PHL 1811 continuum compared to those from gas illuminated by 
the K97 continuum. The ratio is smaller for the high-ionization lines, as might be expected; 
those ions are few in the gas illuminated by the PHL 1811 continuum, so the high-ionization 
line emission is weak. However, for the low-ionization lines, the situation is much different, 
and in fact lines from the same ion exhibit different ratios. Some lines are stronger in 
PHL 1811; for example, Si II A1194 is predicted to be 3.3 times larger for PHL 1811. Other 
lines from the same ion are weaker in PHL 1811; for example. Si II A1814 is predicted to 
be 1.7 times weaker. This difference is a result of the different excitation mechanisms for 
the hues. Mg II A2800 and Si II A1814 are transitions from the lowest permitted levels; 
they are coUisionally excited and their emission strengths depend on the temperature of 
the gas. Since the temperature is lower in gas illuminated by the PHL 1811 continuum, 
the emission from these lines is weaker. In fact, it can be shown that the ratios of the 
emission rates produced by the PHL 1811 continuum and the K97 continuum as a function 
of depth correspond fairly closely to the ratios one would expect based on the difference 
in temperatures through the Boltzmann excitation equation. This is demonstrated by the 
close correspondence between the solid grey hne and the dashed grey hne for Si II A1814 
in Fig. 21 for large depths where Si+ dominates. In contrast, high-excitation Si II, such as 
Si II A1194, is predominately excited by pumping by the continuum. As noted in §4.1,1, 
the UV continuum will be relatively strong in PHL 1811 for the same ionization parameter 
and therefore, lines that are excited by continuum pumping can be expected to be relatively 
strong. This is demonstrated by the great discrepancy in the ratio based on the Cloudy 
emissivity output (solid black line) and the ratio based on the Boltzmann equation (dashed 
black line) in Fig. 21. 

Fig. 20 shows that several bands of Fe II behave differently for the PHL 1811 continuum 
compared with the K97 continuum. We discuss the Fc II further in §4.4. In addition, we 
find that the Fe"*"^ ion column is elevated in PHL 1811, similar to the other intermcdiate- 
ionization hne columns, suggesting that Fe III should be strong. We cannot check this 
suggestion, however, because Cloudy does not model UV Fe III lines, such as Fe III UV 34, 
near 1900A, or Fe III UV 68, near 1950A. 

It is important to note that common UV semiforbidden lines (C III] A1909, Si III] A1892, 
N III] A1750, C II] A2325, N II]A2141) are suppressed for the PHL 1811 continuum. This is a 
consequence of two factors. For the twice- ionized ions, with the exception of Si"^^, the column 
of the ion is lower in PHL 1811. In addition, these lines are all coUisionally excited, arising 
from ions that have two electrons in their P-shell. Thus, since the gas temperature is lower 
in PHL 1811, these lines are weaker. This result has a consequence for plasma diagnostics. 
Semiforbidden lines are generally used as density indicators, since they have critical densities 
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in the range thought typical of broad-hne region gas (e.g., for the hnes above, the critical 



densities range from 3.16 x lO'^cm"'^ for C II A2325 to 3.12 x 10^^ cm"'^ for Si III A1892; 



Hamann et al. 1120021 ). This means that the absence of these lines in objects with very soft 



continua like PHL 1811 does not require the presence of high-density gas. 



4-1 .4- Exploring the logf/, \ogn Parameter Space 

The discussion above considers only a single ionization parameter and density; next 
we examine the properties of the g function of these two parameters. We perform 

simulations for the PHL 1811 continuum and for the K97 continuum for ionization parameter 
—3.5 < logf/ < 0.5 and for density 8.5 < logn < 12.5. In each case we adjust the column 
density so that logiVj^ — logU = 26; this ensures that the gas is very optically thick to the 
continuum, and probes approximately the same thickness relative to the hydrogen ionization 
front. 

We first examine the column density of various ions as a function of ionization parameter 
and density. Since the total column density of the gas varies with ionization parameter, we 
look at the column density relative to the total column; this is equivalent to log(A^jon) — log U. 
Generally speaking, the behavior of the ionic column density for low and intermediate ion- 
ization metal ions relative to the total column density was divided into two groups according 
to ionization potential. For the singly- ionized ions C"*" and N"*", the relative column was 
principally density independent, and decreased as ionization parameter increased, as the 
ionization state of the gas shifted to a higher ionization state. O"*" looked like H+, most 
probably because the strong influence of the Bowen effect. In contrast, the ions Mg+, A1+, 
Si"*" and Fe"*", were almost independent of ionization parameter and density. In both cases, 
the ionic structure was nearly the same for the PHL 1811 continuum and the K97 contin- 
uum; this is not surprising because they should both be roughly equally efficient in producing 
low-ionization ions. 

For the twice-ionized ions, there are significant differences between the gas illuminated 
by the PHL 1811 continuum and the K97 continuum. First, there were differences in the 
total column densities. For the doubly- ionized ions with higher ionization potential (C^^, 
N"*"^, and O^^), the columns were higher for the K97 continuum. In contrast, for the doubly- 
ionized ions with lower ionization potential (Al'*'^, Si"*"^, and Fe"*"^), the columns were higher 
for the PHL 1811 continuum. This makes sense, and it is something that we noted previously 
in CLE: for softer continua, the cooling in the region of intermediate ionization shifts to ions 
that have lower ionization potential. For these ions, the dependence on ionization parameter 
and density were similar: there is an ionization parameter at which the column is maximized. 
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and at that ionization parameter, the column decreases as a function of density. This happens 
because at higher densities, recombination is stronger, and the overall ionization state of the 
gas is a little lower. At lower densities for the same ionization parameter, the ionization 
shifted to the next higher state. But for the higher-ionization doubly-ionized ions, that 
structure was apparent only for the K97 continuum; for the PHL 1811 continuum the column 
increases with ionization parameter almost independent of density. 

More highly ionized ions (C+'^, N+^, He^^) were significantly deficient over the entire 
parameter space for the PHL 1811 continuum, as expected. However, these ions were present 
at highest ionization parameters. 

We next look at equivalent widths of lines that are common in AGN. Note that by 
discussing the equivalent width rather than the line flux, we account for the fact that the 
optical and UV continuum against which the lines are measured is brighter by a factor of ~ 
3.5 in PHL 1811 for the same ionizing flux. We first looked at lines that had equivalent widths 
using the PHL 1811 continuum that were predicted to be larger than or comparable to those 
obtained using the K97 continuum, both at our fiducial point (logC/ = —1.5 and logn = 11), 
and over the parameter space. The strongest hues were the UV Fe II and the high-ionization 
Si II. These are both among the strongest lines in the observed spectrum, a fact that suggests 
that we are on the right track. But not all low-ionization lines were strong in PHL 1811; 
notably weak ones were Mg II and Si II A1814. The difference again is the excitation 
mechanism: the high-excitation silicon and iron lines arc pumped by the continuum, while 
the magnesium and low-excitation silicon are collisionally-cxcitcd. In addition, as expected, 
the high ionization lines such as C IV were among the weakest in PHL 1811 compared with 
K97 continuum. We illustrate this in Fig. 22. At the fiducial point, the equivalent width of 
C IV is a factor of more than 100 times lower for the PHL 1811 continuum compared with 
the K97 continuum. 

The soft continuum also affects the semiforbidden lines in interesting ways. Fig. 23 
shows the equivalent width contours of Si III] A1892 for the PHL 1811 continuum and the 
K97 continuum. This emission line has a critical density of 3 x 10^^ cm~^, so the decline 
toward higher density at about this value for the K97 continuum is expected. But the 
PHL 1811 continuum simulations show a different behavior; the line begins to decrease at 
much lower densities, around 10^°cm~^. Again, this means that the usual density indicators 
fail for soft continua, and the absence of semiforbidden lines does not require the presence 
of high-density gas. 

The weakness of the coUisionally-excited metal lines in PHL 1811 has important con- 
sequences for the physical conditions in the gas. If the gas cannot cool by emission of 
coUisionally-excited lines, and yet there are plenty of photoionizing photons, the energy goes 
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into the production of hydrogen and other continua. This is wh at we see in PHL 1811 



The s ituation is very similar to high-density gas (discussed by, e.g.. lRees. Netzer fc Ferland 



19891 ). At high density, colhsional deexcitation of metal ions prevents the gas from cooling 
by coUisional-excitation of metal ions, and the high-density gas also becomes dominated by 
continua. In a sense, both types of gas are "cooling challenged". They also have similar 
observational properties, in particular, the weakness of the semiforbidden lines. 

To summarize this section, we find that properties of metal ions and their emission in 
gas illuminated by the soft PHL 1811 continuum is different from gas illuminated by the 
K97 continuum in the following ways. The column densities and line emission from high- 
ionization ions are lower in gas illuminated by the soft SED. The column density of the 
intermediate-ionization ions is larger, but the emission from those ions is lower if the line 
is collisionally dominated, or higher, if the line is pumped by the continuum. Thus, the 
strongest lines should be Fe II and high-excitation Si II, as observed, while high-ionization 
lines and semiforbidden lines are weak. As a consequence, semiforbidden lines are not good 
density diagnostics in gas illuminated by the soft SED. Line emission from gas illuminated 
by a soft SED has many similarities with high-density gas, and in principle may be mistaken 
for high-density gas, as a consequence of inefficient cooling by line emission in both cases. 



4.2. A Locally Optimally-emitting Cloud Model 

In the previous section, we showed that the soft SED from PHL 1811 results in low 
emission line flux and low equivalent widths compared to the more typical SED used by K97. 
We first showed this for the fiducial parameters (}ogU = —1.5, logn = 11, logNff = 24.5) 
and then for a range of ionization parameters and densities, with the column density adjusted 
so that log A^i^ + logf/ = 26. As a final investigation into this point, we compute a Locally 
Optimally-emitting Cloud (LOC) model, to ensure that our result is robust against averaging. 



The LOC model was first presented by [Baldwin et al. I (119951 ). Motivated by the fact 



that it is unrealistic to expect that the net emission from the broad-line region could be from 
gas characterized by a single ionization parameter, density, and column density, they postu- 
lated that the observed emission should be an average of emission from gas characterized by 
a range of properties. In the standard LOC model, the radial distribution (equivalent to pho- 
toionizing fiux) and density distribution of the line-emitting gas are characterized by power 
laws, with the radial and density indices equal to —1, and the column density set equal to 
10^^ cm~'^. These parameters are rather arbitrary, but they seem to give relatively good fits to 



most q uasars, at least for the high-ionization lines. For example, ICasebeer. Leighly &: Baron 



( 120061 ) computed an LOC model for RE 1034+39, an object known for its hard SED, and 
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found that the ratio of the prominent hnes with Lya matched the observed ones well; the 
exceptions were C IV and Mg II, which were both predicted to be brighter with respect to 
Lya than observed. 

To compute the LOG model, we run Cloudy models for the following parameters: pho- 
ton flux 17 < log$ < 24, with A$ = 0.125; density 7 < logn < 14 with An = 0.125; 
col umn density Nu = 1 0^^ cm~^ uniformly. This is the range of parameters computed 



by iKorista et al. I (119971). and extends the range originally used by [Baldwin et al. I (119951 ) 



and ICasebeer. Leighly fc Baron I (120061 ) to higher photon fluxes and lower densities. We 
compute the LOG assuming the power law indices of —1 for the radial and density dis- 
tributions; again, these are the standard choices. The outputs of the LOG program are 
the ratios of the line fluxes to the continuum at 12 15 A. We use the input continua to con- 
vert these values to equivalent widths. We choose a covering fraction of 0.25 because that 
seems to give values for the K97 continuum that correspond relativel y well with equ i valent 
widths of genera 
(LBQS quasars), 



ly bright lines from compo site spectra compi l ed by iFrancis et al. I (jl99ll ) 



Zheng et al. I (Il997f ) fi/gr). iBrotherton et al. I (120011 ) {FIRST smvey), and 



Vanden Berk et al. I (l200ll ) (SDSS quasars). The results are given in Table 3. 



Examination of the table shows that the LOG model for the K97 continuum gives a 
reasonable match to composite spectra equivalent widths overall, within a factor of two 
for most of the lines. For example, it predicts a G IV equivalent width of 47A, while the 
four composite spectra yield measured values of G IV of 37, 59, 33, and 24A. This good 
correspondence is expected, since the K97 continuum is intended to be representative of the 
continuum of an average AGN, and the LOG model was developed originally to explain the 
emission lines of composite spectra. 

The equivalent widths of the lines observed in PHL 1811 are much smaller than observed 
in the composite spectra. For example, the observed equivalent width of G IV in PHL 1811 
is only 6.6A, while the composite spectra yield 37, 59, 33, and 24A. Gorrespondingly, the 
LOG model for the PHL 1811 continuum yields consistently much smaller equivalent widths 
for all lines, in some cases by large factors. For example, the equivalent width of G IV for 
the PHL 1811 LOG model is 0.9A, compared with 47A for the K97 LOG model. In fact, the 
PHL 1811 LOG equivalent widths are generally smaller than the observed equivalent widths 
by large factors; the exception is Mg II in which the PHL 1811 LOG model gives an equivalent 
width of 15 A, while the observed equivalent width is 13 A. The largest discrepancies are Lya 
and N V. This may be evidence for a separate optically-thin component that is not well 
described by the standard LOG parameters. We note that difficulties in using the standard 
LOG distributions have been r eported for SPSS J1 54651. 75-1-525313.1, an object that also 
has strong low-ionization lines (jPhanda et al. II2007I ). 
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Despite disagreement between the detailed predictions of the LOG models and the data, 
we have shown that the equivalent widths of the lines in gas illuminated by the PHL 1811 
continuum remain low even when we average over a large range of ionization parameters 
and densities. In the spirit of the locally optimally emitting cloud concept, we note that the 
extremely soft continuum of PHL 1811 allows almost no cloud to emit lines very optimally, 
since gas illuminated by the soft SED has trouble producing line emission. However, even 
though the continuum may be relatively strong in the gas illuminated by the PHL 1811 
continuum (as discussed in §4.1.2), we do not expect to see strong continuum emission. The 
ratio of the Balmer continuum flux to the continuum at 3645A (for the 0.25 covering fraction) 
is 310A for the K97 continuum, but only 94A for the PHL 1811. 



4.3. Fe II and Fe III Emission 

In §3.2.1 we showed that the I Zw 1 Fe II template does not describe the Fe II and Fe III 
pseudocontinuum in PHL 1811. Thus, it is clear that understanding the Fe II and Fe III 
spectrum of PHL 1811 is important for understanding the physics of the line-emitting gas. 
However, the tools that we can use are limited. First of all. Cloudy does not include any near- 
UV Fe III lines. Second, the Fe"*" ion in Cloudy has only 371 levels and the highest one has an 
energy of just 11.6 eV; in §3.2.1 we discussed why we believe that significant emission arises 
from the levels between 13 and 14 eV. In addition. Cloudy uses an approximation for the 
radiative transfer; this might be expected to be a serious model weakness especially for the 
emission from the partially-ionized zone where radiative transfer can influence the thermal 
equilibrium. Furthermore, it has been shown that this Fe II model does not explain the 



spectra of AGN well at all (e.g.. [Baldwin et al. II2004I ). These are clearly serious inadequacies 



for understanding the Fe II and Fe III in this object; nevertheless, we can learn some very 
useful things using the Cloudy Fe II model. 

We start with our baseline parameters (logf/ = —1.5, logn = 11, and logNu = 24.5), 
and run Cloudy models for the range of spectral energy distributions discussed in §4.1.1. We 
first examine the Fe II flux in each of several bands: 1000-2000A (far UV), 2000-3000A (UV), 
4000-6000A (optical), and 7800 A-3 microns (infrared). As can be seen in Fig. 24, in each of 
these bands, the flux from the gas illuminated by the PHL 1811 is second- greatest, following 
and comparable to gas illuminated by the kT = 10 eV continuum. In order to conveniently 
plot all the results in one panel, we show the ratio of the results with respect to the one from 
PHL 1811. The Fe II flux varies with SED differently in the different bands. In the far-UV 
band, the flux decreases smoothly as the SED becomes harder. For the other three bands, 
the flux flrst decreases as the flux becomes harder, but then recovers for hardest SEDs. This 
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effect is most pronounced for the optical band, in which the flux for the kT = 320 eV SED 
is almost the same as that from the PHL 1811 SED. 



As discussed in §3.2.1 and in previous papers (e.g. jDietrich et al. 112002 ), the Fe II/Mg II 



ratio i s used as an abundance diagnostic in quasars, although as discussed in 



Leighly fc Moore 



(120061 ). an observed range in this ratio over a limited range in redshift observed in SDSS 
quasars makes this usage somewhat questionable without a better understanding of Fe II 
excitation. In the second panel of Fig. 24 we show the ratio of the Fe II model emission 
between 2200 and 3050A to Mg II. This ratio is rather constant for harder SEDs, but is very 
high for PHL 1811, a factor of 5.2 higher than that from K97. This very high ratio is a 
consequence of both the higher flux in Fe II (a factor of 2.1 higher in PHL 1811 compared 
with K97) and a lower flux in Mg II (a factor of 2.4 lower in PHL 1811). We have already 
discussed in §4.1.3 the fact that Mg II is a collisionally-excited line and is thus suppressed 
in PHL 1811 because of the low electron temperature. 

As noted in §3.2.1, PHL 1811 does not have a particularly high ratio of optical Fe II to 
H/? compared with a previously-analyzed sample of NLSls; the optical Fe II emission is not 
as spectacularly prominent in PHL 1811 as it is in some AGN. Indeed, the model predicts 
only a factor of 1.7 higher flux in the optical band in PHL 1811 compared with K97, in 
contrast to the factor of 2.2 higher flux in the UV band. Thus, the observed relatively high 
UV-to-optical flux ratio is qualitatively predicted. 



Baldwin et al. I (120041 ) discuss the properties of the 371-level Fe II atom currently imple- 
mented in Cloudy. The model does not fit the observed spectra from AGN at all well in the 
near UV. Specifically, strong emission from resonance lines and low-lying levels near 2400A 
and 2600A are predicted to produced strong "spikes" that are not seen in AGN spectra. 



Baldwin et al. I (120041 ) parameterize this problem using a "spike/gap" flux ratio parameter. 



The parameter is defined as the flux in the 2312-2328 and 2565-2665A ba nds (the spikes 



divide d by the flux in the 2462-2530A band (the region between the spikes). [Baldwin et al. 



( 120041 ) note that the observed ratio is typically near 0.7 in AGN. We plot the spike/gap 
parameter in the bottom panel of Fig. 24. We find that this parameter is minimized for the 
PHL 1811 continuum, so although it is still much larger than observed, the trend is in the 
right direction. 

We can understand these results using the supplemental output from Cloudy which 
includes the departure coefficients and populations of a few selected levels as a function of 
depth, the column density of the 371 levels, and the emissivity of various spectral bins as a 
function of depth. The departure coefficients show that the level populations are consistent 
with those predicted by the gas temperature for upper levels less ~ 7eV. Referring to the 
diagram shown in Fig. 6, these are the upper levels of the resonance lines and low- lying 
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transitions. For higher levels, the departure coefficient for PHL 1811 is consistently higher 
than those from the harder SEDs, and can be very large in the partially-ionized zone. Fig. 
25 shows one of the more extreme cases. High levels are apparently being selectively excited 
in PHL 1811, and since we see an inverse dependence with the hardness of the SED, it is 
most probable that the incident SED is responsible for the pumping. 

The departure coefficients are high in gas illuminated by the PHL 1811 continuum, 
but the temperature in the partially-ionized zone is low, so we need to also examine the 
populations in the levels as well. For levels with energies less than ~ 7eV, the populations 
in PHL 1811 were lowest, and the populations increased as the SED became harder. This is 
expected from the departure coefficient behavior, as the temperature in the partially-ionized 
zone is highest for the hardest SEDs. But for higher levels, the populations reverse, finally 
becoming much larger for PHL 1811 than for the other SEDS. An extreme example is seen in 
Fig. 25. Thus the departure coefficient is high for PHL 1811 both because the temperature 
is low, and because the density of highly-excited ions is high. 

Examination of the populations and departure coefficient shown in Fig. 25 reveals an- 
other important feature of gas illuminated with the PHL 1811 continuum. As will be dis- 
cussed in §4.4, the transition from He"^ to neutral He coincides with the hydrogen ionization 
front in gas illuminated by a typical AGN continuum. At the same time, Fe"*" is almost ex- 
clusively found in the partially-ionized zone. However, for soft continua, the Hc+ to neutral 
He transition occurs at a shallower depth than the hydrogen ionization front. Thus, there 
exists a region in which hydrogen is fully ionized and He is neutral. Fe+ is present in this 
region where hydrogen is fully ionized. This is the origin of the shelf between 8 x 10^° and 
1 X 10^^ cm seen in Fig. 25. These Fe"*" can contribute significantly to the emission. While 
this is not a very extensive region compared with the whole partially-ionized zone, it is worth 
remembering that the UV Fe II is highly trapped, and the bulk of the emission comes from 
the vicinity of the hydrogen- ionization front. The density of Fe"*" in the H^-neutral He zone 
is small, and thus they may be less trapped. Examination of the emissivity times depth 
confirms that a substantial fraction of the emission comes from this region (Fig. 26). 

We are now in the position to understand the results shown in Fig. 24. High-energy 
lines in the far UV are purely pumped; thus they are strongest for the PHL 1811 continuum 
and other soft SEDs, and are weakest for the hardest SEDs. The UV lines have contributions 
from pumping in the higher-energy transitions, and contributions from collisional excitation 
for the resonance and low-lying transitions. Thus, they are strong for the soft SEDs, where 
pumping dominates, and strong for the hardest SEDs that produce the highest temperatures, 
allowing collisional excitation to become equally important. Optical emission comes from 
transitions from the UV levels near 5 eV down to the metastable levels near 3eV; these 
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low-lying levels are predominately coUisionally excited. Recall that the optical emission 
is not trapped like the UV emission, so the hotter partially-ionized zone produced by gas 
illuminated by the hard SED produces plenty of optical Fe II. Finally, the infrared emission 
is produce by Lya pumping; it is thus approximately proportional to the diffuse continuum, 
which is shown in §4.1.2 to be strongest for the soft SEDs, becoming weaker for intermediate 
SEDs, and finally increasing again for the hardest SEDs. 

Overall, the strong UV Fe II emission seems to come from continuum and Lya pumping 
and a significant column of Fe+ in the H+-neutral He zone. All levels of the 371-level atom 
arc more highly populated in the gas illuminated by the PHL 1811 SED. At the same time, 
because of the lower temperature, coUisional excitation is less important, producing less op- 
tical emission relative to UV in the partially-ionized zone. Mg II is also weak, as discussed 
in §4.1.3, and therefore the observed relatively high UV Fe Il/Mg II and UV/optical ratios 
are explained. Finally, the "Spike/Gap" result is explained. The spikes are predominately 
coUisionally excited, and are thus very strong for the hardest SEDs, while the gap is predom- 
inately emission from high levels, and are thus stronger for the softer SEDs. Thus the softer 
SEDs produce a lower "Spike/Gap" ratio, although for the fiducial choice of parameters is 
it still much larger than typically observed in AGN. 

So far we have only examined the Fe II and Fe III emission for a single value of the 
ionization parameter and density. Next, we investigate the effect of varying the ionization 
parameter and density. As was done in §4.1.3, we compute the Cloudy models for the 
PHL 1811 continuum and the K97 continuum for —3.5 < logf/ < 0.5, and 8.5 < logn < 12.5. 
The column density was adjusted so that logiVj^ -|- logU = 26. First, as we noted above. 
Cloudy does not output any Fe III lines; however, we find that the column density of Fe"^^ 
ions is consistently larger in the gas illuminated by the PHL 1811 continuum by a factor of 
~ 2.5, with the largest column densities produced at the highest densities and logC/ = —0.5. 
This is evidence that the Fe III flux should be greater from gas illuminated by the PHL 1811 
continuum. 

We next examine the dependence of the Fe II properties considered in Fig. 24 as a 
function of ionization parameter and density for gas illuminated by the PHL 1811 and 
the K97 continua. For the Fe II flux, we see the same trends as seen in Fig. 24: the 
flux is uniformly larger for the gas illuminated by the PHL 1811 continuum. The largest 
PHL 1811/K97 flux ratio is seen in the FUV band, ranging from ~ 1 to 4.2 with a mean of 
2.9, with the largest ratios observed at the highest ionization parameter and density. The 
same pattern is seen in the UV band, ranging from ~ 1 to 2.6, again with the largest values 
seen at highest ionization parameters and densities. The optical flux ratio ranges from ~ 1 
to 1.65, with a mean of 1.3, and the IR flux ratio ranges from ~ 1 to 2.7 with a mean of 
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1.8. The maximum ratios for the optical and IR ratios occurs however at \ogU = —1.5 and 
logn = 11. However, perhaps a more relevant parameter is the equivalent width, since for 
any photon flux, the PHL 1811 continuum will be brighter than the K97 (Fig. 13). The 
equivalent widths in the far UV and IR bands are almost the same, and they are a factor of 
~ 2 lower in the UV and optical. 

We next look at the behavior of the ratio of the UV Fe II flux in the 2200-3050A band 
to the Mg II flux as a function of ionization parameter and density (Fig. 27). Generally 
speaking, this ratio is larger for gas illuminated by the PHL 1811 continuum compared with 
the K97 continuum by a factor of ~ 3 or more. As can be seen in Fig. 27, the range in 
the ratio is large with a strong dependence on ionization parameter; however, the mean 
obtained from the log of the ratio is ~ 2.8 for the K97 continuum, which is not far from 
the value measured from composi te spectra from the SDSS (3.9 or 5.3, depending on Fe II 



template, iLeighly fc Moore II2004I ). The mean ratio for the PHL 1811 continuum is much 
higher, around 7.9. The ratio is larger for the PHL 1811 continuum, both because the UV 
Fe II flux is higher by a factor of ~ 1.8, and because the Mg II flux is lower by a factor of 
~ 1.5. 

Finally, we investigate the dependence of the "spike/gap" ratio as a function of ionization 
parameter and density (Fig. 28). This ratio is consistently smaller for gas illuminated by 
the the PHL 1811 continuum; the mean for the PHL 1811 continuum is 3.5, and the mean 
for the K97 continuum is 4.6. Thus the excitation differences discussed are maintained as 
the ionization parameter and density is varied. We note, however, that even though the 
spike/gap ratio is lower for PHL 1811, the min imum is 2, which is still much larger than the 



observed value of ~ 0.7 (IBaldwin et al. II2004I ). But as discussed above, the relatively larger 



populations at high levels in the gas illuminated by the PHL 1811 is an indication that a 
more complete atom would produce emission from higher levels, plausibly filling in the gap 
further. 

To summarize this section, we note that the prominence of the Fe II and Fe III in 
PHL 1811, and the unusual shape of the pseudocontinuum mean that understanding the 
emission of these lines is important for understanding PHL 1811 in general. Although the 
tools at our disposal are limited ( Cloudy does not predict U V Fe III emission lines, does not 
do radiative transfer properly, and its Fe"*" model atom only goes up to 11.3 eV), we show 
that the soft SED can qualitatively explain features of the Fe II and Fe III emission. A 
higher predicted column density of Fe"*"^ predicts stronger Fe III, as observed. Fe II emission 
is predicted to be generally stronger in gas illuminated by the soft SED, as a consequence 
of the continuum pumping and interaction with the strong diffuse continuum. This trend is 
important in the near UV, resulting in a high Fe II/Mg II ratio, as observed. The optical is 
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less affected, resulting in a high UV-to- optical Fe II rat i o, as observed. High excitation and 
departure coefficients, as well as a low [Baldwin et al. I (120041 ) "spike" to "gap" ratio, show 
that emission from high levels is stronger in gas illuminated by the soft SED, suggesting that 
the observed strong emission from very high levels may be predicted if a larger model atom 
were available. 



4.4. "Filtering" and He I, Ha, and narrow Lya 

As we show in §3.1, we do not detect He I A5875 in PHL 1181. Assuming that He I A5875 
has the same profile of Ha, we obtain a conservat i ve upp er limit of 0.0062 on the He I/Ha 



ratio in PHL 1811. As noted before, I Crenshaw I (119861 ) tab ulated the fl u xes o f He I and 



Ha for 9 Seyfert galaxies; the mean He I/Ha ratio is 0.061. [Thompson I (jl99ll ) tabulated 
He I/Ha ratio for 6 quasars, and the mean ratio is 0.025. Thus, PHL 1811 has abnormally 
weak He I emission. 

Before going any farther, it is worthwhile reflecting that He I and Ha are recombination 
lines and are thus considerably different from the lines that we have discussed already. Also, 
in gas illuminated by a typical AGN SED, the transition from He"*" to neutral He occurs 
at approximately the same depth as the transition from H+ to neutral H - the hydrogen 
ionization front. In a gas illuminated by a soft SED, the continuum runs out of photons that 
can ionize helium before it runs out of photons that can ionize hydrogen. The result is that 
there is a region in which hydrogen is ionized and helium is neutral. He I is produced in this 
region, which means a low He I/Ha ratio is more difficult to produce for soft SEDs than for 
hard ones. 

We first use the Cloudy simulations using the CLB continua to explore whether we can 
produce a low He I/Ha ratio, especially for soft continua. In this set of simulations we exam- 
ine the ratio of the simulated flux to the observed flux evaluated at Earth. Representative 
results as a function of kT^ut, ionization parameter, density and depth D, where depth is 
measured in terms of the depth of the hydrogen ionization front are shown in Fig. 29. We 
find that Ha is almost independent of kTcut and increases with column density. Ha decreases 
at high density at a given ionization parameter because of thermalization. We find that we 
can produce the Ha emission that we see for kTcut = 10 and D = 10 when log(f/) = —2 and 
log(n) = 11.5, although we note that we have assumed 100% covering fraction, and including 
a more realistic covering fraction would reduce the model flux. In this region of parameter 
space, and every other where we can match the Ha emission, the predicted He I emission 
is far too strong. He I decreases with density; however. Ha decreases more quickly, and at 
very high densities and ionization parameters. Ha is predicted to be too weak. As discussed 
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above, we see that at a given depth the He I flux increases as the SED becomes softer. 

It is possi ble to ob t ain a low He I/Ha ratio using the concept of "filtering" that was first 



introduced in iLeighly I (120041 ) . In that paper, it was shown that the intermediate- and low- 
ionization lines of two NLSls could be better explained if the continuum were transmitted 
through the wind creating the blueshifted high-ionization lines observed fro m these objects 



before it illuminates the intermediate- and low-ionization line-emitting gas (see lLeighly II2004J . 
for details). 

PHL 1811 has blueshifted high-ionization lines, and perhaps filtering can help with the 
He I/Ha problem. Specifically, the continuum may be transmitted through a wind so that 
the emerging spectrum lacks photons in the He I continuum and yet the column density is 
such that the hydrogen continuum remains relatively strong. In fact, it is easier to produce 
such a filtered continuum in gas illuminated by a soft SED, because of the presence of the 
region of gas in which hydrogen is ionized and helium is neutral. The presence of this region 
means that the continuum will run out of photons in the He I continuum before it runs out 
of photons in the hydrogen continuum, yielding the transmitted continuum that we need. 

We test this scenario using Cloudy. We first show that a continuum with the necessary 
features can be produced. We use the PHL 1811 continuum shown in Fig. 12. The high 



N V/H e II ratio implies that the metallicity in PHL 1811 may be enhanced; as in iLeighly 



( 120041 ). we assume metals augmented by a factor of five, and nitrogen enhanced by a factor 
of ten. We first locate the column density at which the ratio of the He+ column to the 
H^ column is the largest for a range of densities 7.875 < log(?2) < 12.75 and ionization 
parameters —3.0 < log(t/) < —1.125. We then run Cloudy models using these column 
densities, and output the transmitted continua. We identify SEDs with weak He"*" continua 
by taking the ratio of the transmitted plus diffuse emission at 23.9 eV, just below the He I 
edge at 24.6 eV, and at 28.8 eV, just below the N II edge at 29.6 eV, which seems to be 
the most prominent feature toward higher energies from the He I edge. We find a very large 
range in this ratio, from 2.6 to larger than 10^. In fact, there is a plateau in the ionization 
parameter /density plane such that the ratio high for densities less than about 10^^'^, and 
for ionization parameters higher than —2.0 (Fig. 30). We choose a representative continuum 
with log{U) = —1.5 and log(r;,) = 11.5, where the column density at which He"'"/H"^ is 
maximized is lO^^'^cm"^. We note that any continuum from the afore-mentioned plateau 
should produced the same result. The continuum is shown in Fig. 30. 

The gas that the continuum is transmitted through will produce line emission. This 
gas may be th e wind that emits the blueshifted lines, as was assumed and required in 



Leighly I (120041 ). The Cloudy simulations show that the strongest predicted line relative to 



the observed emission lines is Lya. A covering fraction of 10% ensures that it is predicted 
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to be no stronger than observed. Also, if the gas filtering the continuum is a wind, the 
line emission may be broadened and blueshifted, and would be more difficult to distinguish 
from the continuum than a narrow line. As shown in Fig. 30, much of the energy of the 
absorbed continuum is emitted as bremsstrahlung in the infrared. This emission will join 
the transmitted continuum in illuminating the low-ionization line-emitti ng gas, where it may 



enhance the low-ionization line emission through free-free absorption (IFerland fc Persson 



19891). 



Finally, we run Cloudy simulations to determine whether the predicted He I is reduced 
as expected; the results for He I are shown in Fig. 31, and the results for Ha are shown in 
Fig. 32. Each figure shows the predicted emission compared with observed as a function of 
ionization parameter and density for the unfiltered and filtered continuum shown in Fig. 30. 
We assumed a covering fraction of 10%, noting that the covering fraction must be equal to or 
smaller than the one required by the line emission in the filtering gas, above. These figures 
show that Ha is predicted to be nearly the same for the direct PHL 1811 continuum and 
for the transmitted continuum shown in Fig. 30. This is not surprising because the Lyman 
continuum, responsible for photoionizing the hydrogen, is not much different for the direct 
and transmitted continua. Including the covering fraction means that we produce sufficient 
Ha for moderate-to-low ionization parameters, high column densities and moderate densities, 
although we are a factor of two lower than the observed. In contrast, the He I is much weaker 
for the transmitted continuum, and is now in the vicinity of the measured upper limit for 
the same parameter choices. 

To summarize this section, we used Cloudy models to try to understand the very low 
He I/Ha ratio observed from PHL 1811. We found that a soft spectral energy distribution 
exacerbates the problem; He I should be stronger in objects with soft SEDs because of the 
large H+/neutral He zone. We postulated that transmitting the continuum through a wind 
may yield a "filtered" continuum with a large He I edge, noting that a similar scenario 
was successful in e xplaining the in termediate- and low-ionization lines in IRAS 13224—380 



and IH 0707—495 (iLeighly II2004J ). To test this idea, we showed that for column densities 
chosen to maximize the He+ column to the H+ column ratio, a large area of log(f/)-log(n) 
parameter space would produce transmitted spectral energy distributions with the He I 
continuum suppressed by a factor of up to 10^. That we can create such a continuum 
is a direct consequence of the presence of a region in gas illuminated by a soft spectral 
energy distribution where fully ionized hydrogen and He"*" coexist; it means that the column 
density of the filtering gas can be chosen to lie in this region, producing a large He I edge 
in the transmitted continuum. We show also that transmission through such gas would 
not be predicted to produce strong line emission that is not observed if we assume a 10% 
covering fraction. Rather, for a representative case, we see that most of the energy goes 
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to bremsstrahlung emission in the infrared. We illuminate gas with the filtered continuum, 
and find that while Ha remains very similar to the unfiltered case, He I is significantly 
suppressed. Thus, we have demonstrated again that some line emission in AGN can be 
explained by illumination by continua that have been transmitted through a wind. 



4.5. Is PHL 1811 Unique? 



PHL 1811, being intrinsically X-ray weak, and lacking strong emission lines, is a very 
unusual AGN. Are there any others like it? We do not know of any exactly like it, but there 
are a few that are similar. 

As seen in Fig. 9, PHL 1092 and RX J0134. 2-4258 have similar far UV spectra; they 
both also lack strong and prominent Lya and C IV lines. PHL 1811 is also similar to 
these objects in that they all have extreme Eigenvector 1 properties; specifically, they have 
very high Fe II/ [O III] ratios, because they all l a ck or h ave relatively weak [O III] emission 



(PHL 1811: Fig. 2; PHL 1092: e.g., iLeighlv I fll999bf ): RX J0134.2-4258: iGrupe et al. 



( 120001 )). They also have unusual X-ray properties. RX J0134.2— 4258 had one of the softest 
X-ray spectra of any AGN observed in the ROSAT All Sky Survey, which implies a steep 
Uox'-i later ROSAT PSPC and ASCA observations found it to have a more typical X-ray 
spectrum and Uox- It showed rather high amplitude X-ray variability during the ROSAT 
pointed observation and the ASCA observation, although typical for a luminous NLSl. Th e 
ROSAT spectrum of PHL 1092 revealed a very steep spectrum (IForster fc Halpern Ill996l ). 
similar to RX J0134.2— 4258. It also showed rat her high amplitude shor t and long time 



scale X-ray variabili t y for such a lumino us quasar (IForster &: Halpern Ill996l : iLeighly 



1999a 



Brandt et al. I Il999l : iGallo et al. I l2004j ). There may be other objects that are similar for 



which we have less information, such as RX J2219.9— 5941, a luminous NLSl which shows 



high amplitude X-ray variability (iGrupe et al. I l200ll . |2004| ) 



The common properties among these three objects are their unusual UV spectrum, their 
hig h luminosity for NLSls (e. g., PHL 1811 is more luminous than all but 2 of the 87 quasars in 
the lBoroson fc Green I (119921 ) subsample of the PG quas ars), their ex treme X-ray properties, 
and their extreme Eigenvector 1 properties. As noted by lLaor I ( l2000l ). high luminosity NLSls 
should have high Eddington ratios. We previ ously speculate d that X-ray weakness in NLSls 
is a characteristic of a high Eddington ratio (ILeighly Il200ll ). XMM-Newton observations of 
four luminous NLSls found that one (RX J1225. 7+2055) was X-ray weak with Uox ~ —2.0 
( iMatsumoto. Leighly &: Kawaguchi II2004I ). Furthermore, a compilation of Uox values from 
NLSls from the literature indicated that X-ray weakness may be more common in luminous 
NLSls. In Paper I, we speculated that a high Eddington ratio leads to a weak corona; the 
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quenching mechanism proposed by iProga I (120051 ) may be applicable 



In the Sloan Digital Sky Survey, some unusual high-redshift QSOs have been discovered 
that have very blue rest-frame UV spectra, but no discernible emission lines. They are 
identified as quasars, and their redshifts are obtained fr om the onset of L ya forest absorption. 
The prototype object is SDSSp J153259. 96— 003944.1 (IFan et al. Ill999l . hereafter referred to 
as SDSS J1532— 0039). Featureless continua are typically found in BL Lac objects. However, 
SDSS J1532— 0039 does not appear to have other characteristic properties of BL Lac objects; 
it was not detected in a deep radio observation and was not strongly op tically polarized. It 
has recently been found to be optically variable (IStalin fc Srianand II2005I ): however, this may 
not be conclusive evidence that it is a BL Lac object, since one luminous NLSl (PHL 1092; 
Gallo et al. 1120041 ) has been observed to be tentatively variable in the UV on short time scales. 
Furthermore, follow-up observation s with Chandra failed to detect the quasar, indicating an 
c^ox < —1.74 (IVignali et al. Il200ll ). We thought that the lack of broad emission lines a nd 
the X-ray weakness suggest a similarity to PHL 1811 (ILeighly. Halpern fc Jenkins 1 120041 ). 



A few other objects appear eq ually puzzling. The z~ 1 quasar PG 1407-1-265 is possibly 



the first known weak-line object (IMcDowell et al. I Il995l ). Two more such high redshift 
obiects were found in SDSS (SDSS J130216. 13+003032.1 and SDSS J144231. 73+011055.3 ; 



Anderson et al. 1 120011 ). The z=0.494 quasar 2QZJ215454.3-305654 flLondish et al. Il2004f ) 



was found in the 2dF/6dF QSO redshift survey. It has a blue continuum, no detected radio 
or X-ray emission, and is not polarized. In the Hamburg/ESO survey, the z=1.8 quasar 
HE 0141—3932 was found to have unusu ally weak Lya emiss ion; It also has blueshifted 
high-ionization lines, similar to PHL 1811 (IReimers et al. II2005I ). 



Because of their steep spectra and high amplitude X-ray variability, it has been pro- 
posed that NLSls are characterized by a high accretion rate. Luminous NLSls should have 
exceptionally high accretion rates compared with their Eddington value. Luminous high 
redshift quasars should have large black holes, and in order to grow so large in the short 
amount of time implied by the standard cosmology, they should be accreting at a rapid rate. 
If the weak X-ray emission and UV line emission is characteristic of a high accretion rate in 
a large black h ole, it may mean that PHL 1811 is the local prototype of the high-2 lineless 
QSOs (see also iLeighly. Halpern fc Jenkins II2004I ). 



The observations, however, do not unequivocally support this argument. First of all, 
not all of the lineless quasars found in the SDSS have weak X-ray emissio n; some have 
normal anm- For example, PG 1407+265 has aox —1.44 (IWilkes et al. Ill999l). whic h is ac- 

' 994h . Also, 



Wilkes et al. 



tually fiat for a quasar of its 2500A luminosity (log(/o) = 32.01 
SDSS J144231. 73+011055.3 has a relatively fiat aox of -1.37 (jSchneiderelal 



20051 ). Fur- 



thermore, it is quite possible that some hne-less objects are radio-quiet BL Lac objects 
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( ICoUinge et al. II2005I ). i.e., a hypothetical class of objects in which the beamed continuum 
overwhelms the lines, but unlike ordinary BL Lacs, there is for some reason weak radio and 
X-ray emission. Thus, it could be that lineless quasars are a heterogeneous population. 

Although there are a number of other objects that are quite similar to PHL 1811, the 
group comprises just a handful, insignificant compared with, for example, the thousands of 
quasars identified in the SDSS. But can we be confident that we would be able to identify 
all of these objects? X-ray emission is considered to be an identifying feature of AGN, and 



it ha s been shown to be useful for identifying AGN in deep fields (e.g., iBrandt fc Hasinger 



20051 ). and is even thought to be a k ey for constrainiii g the AGN contribution to the power 



in ultraluminous /i?74 5" galaxies (e.g.. |Ptak et al. II2003I ). PHL 1811 is a powerful quasar with 



a large black hole, but if it were embedded in a ULIRG it would not be identified as such 
because of its weak X-ray emission. Quasars are also identified by their radio emission, and 
PHL 1811 was rediscovered as a member of the FIRST Bright quasar survey. However, it is 
a weak radio source, and had it been further away, it would not have been found. AGN are 
commonly identified by their blue color in photometric surveys, and PHL 1811 would have 
been easily identified as a PG quasar, if it has been in the region of the sky covered by that 
survey. However, if it were embedded in a dusty starburst host, the continuum would be 
reddened, and it would be much more difficult to find. One way that we might be able to 
identify these objects is through their rest-frame infrared properties. It has been shown that 
AGN can be distinguished from starbursts on the basis of the ^9.5^,,^/ Fm/. m. fiux ratio, which 



measures the temperature of the dust (e.g.. lde Grijp. Lub &: Miley 1119871 ). Otherwise, if we 



speculate that X-ray weakness, and weak line emission is a characteristic of a high accretion 
rate, one might imagine legions of PHL 1811s lying, unidentified, in high redshift galaxies. 



5. Summary 

In this paper, we report the results and analysis of HST STIS optical and UV spec- 
troscopic observations that were coordinated with the Chandra observations discussed in 
Paper I, and contemporaneous ground-based optical observations of the unusually luminous, 
nearby narrow-line quasar PHL 1811. Here, we summarize the primary results of the paper. 

• The optical spectrum indicates that PHL 1811 is a luminous example of a Narrow-line 
Seyfert 1 galaxy, but the emission lines are not typical of any type of AGN. Overall, 
the line emission is weak. There is no evidence for forbidden or semi-forbidden line 
emission. Very low ionization lines are present, including Na ID and Ca II H&K, 
that are not typically observed in AGN. These lines appear to be narrower than H/5. 
The near-UV region of the spectrum is dominated by Fe II and Fe III. The UV Fe II 
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emission is strong with respect to both the optical Fe II and Mg II. We find evidence 
that the UV Fe II and Fe III pseudocontinuum has a different shape in PHL 1811 
compared with other AGN including the prototype NLSl I Zw 1; this appears to be a 
consequence of a higher ionization/excitation state of the iron in PHL 1811. The C IV 
line is remarkably weak, with an equivalent width of only 6.6A, about 5 times weaker 
than C IV in quasar composite spectra; it is also broad and blueshi fted, similar to the 



NLSl s IRAS 13224-3809 and IH 0707-495, indicating an outfiow flLeighlv fc Moore 



20041 ). The region around Lya is strongly blended, but can be successfully deconvolved 
using a template created from the C IV profile into Lya, N V, Si II and C III*. We 
also find that Si III* is clearly present near 1112A. 

In Paper I, we showed that PHL 1811 is intrinsically X-ray weak, and therefore we infer 
that it has an unusually soft spectral energy distribution. Can a soft SED by itself 
produce the observed unusual line emission? We explore this question using Cloudy, 
and find that it can. A soft spectral energy distribution is stronger in the UV for 
a particular photoionizing fiux than a hard SED; this is one factor that reduces the 
observed equivalent widths, but it also enhances lines that are pumped by the contin- 
uum, such as high-excitation Si II. The gas is cooler than gas illuminated by a hard 
SED, but has a higher hydrogen ionization fraction in the partially- ionized zone than 
might be naively expected. The low temperature means that coUisionally-excited line 
emission, including semiforbidden lines and Mg II, is weak. Without the collisionally- 
excited fine emission to cool the gas, the photoionization energy goes to excited states 
of hydrogen and continuum emission, rather than to the electrons, resulting in the cool 
temperature. We conclude that the soft SED alone can qualitatively explain the line 
emission: high-ionization lines are weak because the continuum lacks photons to cre- 
ate the highly-ionized ions, semiforbidden and low-ionization collisionally-excited lines 
are absent because the gas is too cool to excite them, but lines excited by continuum 
pumping are strong because a soft SED continuum is stronger for a given photoionizing 
flux than a typical SED. We also compute an LOG model and show that the above 
results are robust to averaging. 

Fe II and Fe III are very important in the near-UV spectrum. The Cloudy iron ion 
models are too limited to sufficiently model this line emission in detail. Specifically, 
we see strong Fe III in the spectrum, but no UV Fe III lines are output by the model. 
In addition, we infer Fe II emission from highly excited levels near 14 eV, but the 
highest level treated by the model is only 11.3 eV. Finally, Cloudy uses a approximation 
for radiative transfer, and that may be a significant shortcoming in the partially- 
ionized zone. Nevertheless, Cloudy simulations reveal evidence that the Fe II and 
Fe III emission is qualitatively explained by soft SED. The Fe"*"^ columns are larger in 
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gas illuminated by the soft SED, for a range of ionization parameters and densities, 
due to the fact that the H II region is dominated by intermediate-ionization ions, since 
the usual high-ionization ions are not produced by the soft SED. The predicted Fe II 
emission is uniformly higher in gas illuminated by the soft SED. This is a consequence 
of several factors including strong continuum and Lya pumping, as mentioned above, 
as well as the presence of a region of gas in which hydrogen is ionized but helium 
is neutral and where Fe"*" is also present. The excitation of Fe"*" is higher in the gas 
illuminated by the soft SED, and relatively high departure coefficients are seen in the 
highest levels treated by the model, suggesting that emission from even higher levels 
might be strong. Overall, as far as we can tell given the limitations of the model, the 
soft SED is responsible for the unusual Fe II and Fe III observed. 

PHL 1811 has a very low ratio of He I A5876/Hq;. This is not easy to explain using the 
soft SED because a soft SED predicts a larger value of this ratio as a consequence of the 
presence of a region in which hydrogen is ionized but helium is neutral. We observe 
blueshifted high-ionization lines in this object and postulate that if the continuum 
were transmitted through optically-thin gas, the resulting "filtered" continuum may 
be strong in the hydrogen continuu m, necessary to create Ha, but absorbed in the He I 



continuum, with a large He I edge. iLeighly I (120041 ) showed that a "filtered" continuum 



was able to better explain the intermediate- and low-ionization line emission in the 
NLSls IRAS 13224-3809 and IH 0707-495. We locate the regions of parameter 
space for the filtering gas where the required large He I edge will be produced, and 
demonstrate that illuminating gas with the resulting continuum will produce a low 
He I/Ha ratio as observed. 



In this paper, we present a detailed analysis of a single unusual object, PHL 1811. 
However, a number of things that we have learned may carry over to our understanding of 
AGN in general. For example, the PHL 1811 spectrum lacks forbidden and semiforbidden 
lines. Generally, one would interpret weak semiforbidden lines to be a consequence of high 
density. We showed, however, that weak semiforbidden lines are expected from gas illumi- 
nated by a soft SED, which means that semiforbidden lines do not perform well as density 
indicators when the illuminating SED is soft. How many times have weak semiforbidden 
lines been taken to be evidence for high density, when in reality, a soft SED is the culprit? 
The similarity between emission from a high-density gas and gas illuminated by a soft SED 
is not accidental; they are both a consequence of the inefficiency of collisionally-excited line 
emission. In the case of PHL 1811, the soft SED is intrinsic, but it is also possible that the 
soft SED may be the result of filtering. Whether or not a filtered continuum can mimic a 
high density gas will be the topic of a future paper. 
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PHL 1811 is a high-luminosity NLSl, and according to the currently-accepted paradigm, 
these objects should be accreting at the largest rates. Therefore, we may expect that they 
are present at high redshift. However, they may be difficult to identify as AGN, and there 
plausibly be a large number embedded in dusty galaxies. The line emission is weak, and in 
spectra with low signal-to-noise ratio, the optical and UV spectra may appear not to have 
any emission lines at all. The continuum colors are blue, but that would not be seen if 
embedded in a dusty galaxy. The X-ray emission is a hallmark of AGN, but it is weak in 
PHL 1811. The best way to identify similar objects at high redshift may be through their 
infrared properties. 
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A. Appendix 

In §4 we showed that gas illuminated by the soft spectral energy distribution has rather 
unique properties that have not been discussed in the literature in the context of AGN 
emission lines, and these properties produce the unusual line emission in PHL 1811. These 
unique properties include weak high-ionization lines, low temperatures that result in weak 
coUisionally-excited lines, a strong UV continuum for a given value of photoionizing flux 
that results in strong pumped lines, and high level of hydrogen ionization in the partially- 
ionized zone due to the inability of the gas to cool via coUisionally-excited line emission. 
In this appendix we further investigate aspects of the properties of the gas illuminated by 
the soft spectral energy distribution that are not directly related to interpretation of the 
spectra, but are nevertheless interesting and are included for completeness. In particular, 
we are interested in precisely what characteristics of the continuum result in the observed 
properties; is it the weak X-ray or the weak extreme UV, or both? This issue is examined 
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in §6.1. We are also interested in whether we can destroy these characteristic properties by 
including turbulence and changing the metallicity; these issues are discussed in §6.2. 

A.l. Effect of the Spectral Energy Distribution 

The PHL 1811 spectral energy distribution shown in Fig. 12 differs from a typical quasar 
in two respects: it is weaker in the extreme UV, and it is weaker in the X-ray. In fact, we do 
not know the shape of the spectrum in the extreme UV, so in principle it could be brighter 
in that band. In this section, we examine the influence of these two factors by varying the 
extreme UV cutoff and the independently. 

The spectral energy distributions were constructed from the PHL 1811 SED shown in 
Fig. 12, with a few changes. First, the PHL 1811 UV continuum extends to \ogv — 15.53 
before dropping to the hard X-rays. We round this value down to logz/ = 15.5. We then 
create another break in the spectrum, in the extreme UV, and vary the frequency of that 
break between log// = 15.5 and 16.4. The upper limit corresponds to about 100 eV. These 
new points arc extrapolated from logz/ = 15.5 along a power law with a slope (in logvFy) of 
0.5. Then, they are joined to the X-ray, but while the PHL 1811 continuum shown in Fig. 
12 joins the X-ray at logz/ = 16.76, we instead join at 17.1, as that better approximates the 
position of the break in the K97 spectrum. Finally, the X-rays are renormalized to produce 
an between —2.25 (approximately the observed value for weaker X-ray observation of 
PHL 1811) and —1.4 (approximately the observed value for K97). The spectral energy 
distributions are shown in Fig. 33. Cloudy models for these 324 SEDs are run for the fiducial 
parameters (logf/ = —1.5, logn = 11.0, logNn = 24.5), yielding a grid of results as a 
function of aox and the extreme UV cutoff. PHL 1811 is characteristic of steep aox and low 
EUV cutoff (lower left corner of the plots that follow) and K97 is characteristic of fiat aox 
and high EUV cutoff (upper right corner of the plots that follow). 

We first look at the physical characteristics of these spectral energy distributions as a 
function of aox and EUV cutoff. The flux at 2500A is nearly purely dependent on the EUV 
cutoff, with the largest values being found for the lowest cutoff energies; this makes sense 
because the majority of the photoionizing photons have energies just larger than 13.6 eV 
due to the steep dependence of photoionizing flux on energy. Thus, the larger the EUV 
cutoff, the lower the overall normalization for the same photoionizing flux. In contrast, the 
Compton temperature is nearly purely dependent on aox- 

We next look at gas properties as a function of depth and normalized to the depth of 
the hydrogen ionization front. In the H II region, the temperature depends on both agx and 
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the EUV cutoff, with the lowest temperatures found for the steepest aox and lowest-energy 
EUV cutoff, and the highest temperatures for the flattest aox and highest-energy EUV cutoff. 
Then, for lines produced in the H II region, the behavior of the column density of the ion 
producing the line and the lines themselves depends on the ionization potential, primarily, so 
that there is dependence on Uox if the ionization potential is high, and limited to dependence 
on EUV cutoff if the ionization potential is low, roughly speaking. Contours of flux for the 
high-ionization line C IV and the intermediate-ionization lines C III] A1909 and Si III] A1892 
are shown in Fig. 34. 

Turning to the properties of the partially-ionized zone, we examine the temperature, 
the hydrogen ionization fraction, and the ratio of the excitation temperature to the electron 
temperature at n = 2 at the representative depth of ten times the depth of the hydrogen 
ionization front (Fig. 35). The temperature depends primarily on for flatter aox with 
some dependence on EUV cutoff for steeper aox- The lowest temperatures occur, not sur- 
prisingly, for objects with the steepest aox and lowest energy EUV cutoff. The hydrogen 
ionization fraction depends primarily on the EUV cutoff, which shows that photoionization 
from n = 2 by the input continuum is the dominant process. However, at the flattest values 
of ftox, the hydrogen ionization fraction increases somewhat, due to the additional ionization 
by nonthermal electrons from photoionization by soft X-rays. That this depends on aox is 
not surprising, because the break to the X-ray power law, at 100 eV, is in the soft X-ray band 
so that the X-ray power law includes a region of high opacity to photoelectric absorption 
by metals. Finally, the ratio of the excitation temperature to the electron temperature of 
n — 2 hydrogen is a convolution between the temperature and the hydrogen fraction, and 
the highest values occur for the steepest aox and lowest EUV cutoff, similar to the PHL 1811 
continuum. 

The low-ionization lines, excluding Fe II, divide naturally into those that are pumped 
and those that are coUisionally excited. CoUisionally excited lines are weakest for the steepest 
aox and the lowest values of the EUV cutoff. Examples of such lines are Mg II, C II A1335, 
and Si II A 18 14. The pumped lines follow the UV flux and are strongest for the lowest values 
of the EUV cutoff. Examples of such lines are the high-ionization Si II lines including those 
at 1260 and 1194A. Contours for Mg II, Si II A1814, and Si II A1194 are shown in Fig. 34. 
The Fe II hnes are all stronger for low values of the EUV cutoff, but in some cases, they 
increase for flat aox at a particular value of the EUV cutoff. In this way they resemble the 
hydrogen ionization fraction shown in Fig. 35. These may be responding to Lya pumping. 
Examples of Fe II that appears to be pumped are UV 191, the infrared lines between 7300A 
and 3 microns, which are thought arise predominately in cascades after Lya pumping, and 
the "gap" . Example of lines that may be responding to Lya pumping are the "spike" lines 
and other low-lying levels, and the near- UV band between 2000 and 3000A. Fig. 34 shows 
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the results for UV 191, the "spike" and the "gap". 

We also computed the equivalent widths, but since the continuum is much stronger for 
continua with low EUV cutoffs, they are almost all weaker for low EUV cutoffs and steep 
aox- Lines that are purely pumped by the continuum have equivalent widths independent of 
SED shape. 

In this section, we have investigated simulations for only single values of the ionization 
parameter and the density; computing an average or LOG would provide more robust results 
but is clearly beyond the scope of this paper and will be pursued in the future. Nevertheless, 
it seems that we can differentiate between SEDs by the behavior of the high-, intermediate- 
and low-ionization lines. Intermediate-flux high-ionization line emission will be produced for 
either high EUV cutoff or flat ao^^ or both. Weak high-ionization line emission requires both 
a low EUV cutoff and a steep aox- The same is true for intermediate- and collisionly-excited 
low-ionization line emission. Pumped lines are stronger for low EUV flux. 



A. 2. Effect of Turbulence and Metallicity 



In the simulations presented in §4 and in §6.1, we assumed solar metallicity and no 
microturbulence. Both of these factors influence the cooling of the gas, and can therefore 
potentially disrupt the effects of the soft spectral energy distribution discussed previously. 

The most extensive inves tigation of the influenc e of microturbulence on the line emission 
in quasars was performed by iBottorff et al. I (120001 ). They noted that microturbulence has 
two main effects on the spectrum. First, the line optical depths decrease as the line proflle 
is spread over a larger range in wavelength. This causes the lines to be brighter overall as 
they can escape more freely. The second effect is that continuum pumping becomes more 
important, again because the lin e profile has a larger width. As far as the gas properties are 
concerned, iBottorff et al. I (120001 ) indicated that the temperature will be lower for very large 
Vturh as the photons escape more readily and cool the gas better, but higher for intermediate 
values due to continuum pumping of lines. They found that most lines become brighter 
as Vturh is increased, but by the same amount, which means that they would not be good 
diagnostics of the infiuence of turbulence. Some lines, however, b ecome a lot br i ghter , 
including Si II A1260, Si II A1305, C II A1335, and Fe II A1787. Thus, IBottorff et all (l200oh 
suggested that ratios of these lines with respect to intercombination lines and Si II A1814 
would be good diagnostics of turbulence. 



T here have been a num ber of publications that discuss the infiuence of enhanced metal- 
licity. iFerland et al. I (Il996l ) discussed the effects of increasing the metallicity on production 
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of the high-ionization li nes, especially N V, •y y ith th e goal of explaining the high velocity com- 



ponent of Q0207— 398. ISnedden fc Gaskell I (Il999l ) discussed the fact that enhanced cooling 



by higher metallicity gas will result a somewhat higher characteristic ionization parameter, 
allowing a smaller characteristic radius and higher characteristic density, thus alleviating 



discrepancies between photoionization models and reverberation mapping. iHamann et al. 



(120021 ) investigated the behavior of various lines as a function of metallicity, with the aim of 



determining the line ratios th a t are most robustly sensitive to metallicity variations to use as 



diagnostics. Finally, iLeighly I (120041 ) showed that enhanced metallicity is required to explain 
unusual high-ionization line ratios in two NLSls that have very blueshifted high-ionization 
lines. 

In this section, we discuss the effects of the metallicity and turbulence on the properties 
of the gas illuminated by a soft spectral energy distribution. Our goal is qualitative insight, 
so we limit ourselves to the fiducial gas parameters (logf/ = —1.5, logn = 11.0, and logNn = 
24.5). It would be interesting to explore the full parameter space, but that is beyond the 
scope of this paper. 

We first investigate the temperature of the g function of depth normalized to 

the hydrogen ionization front, shown in Fig. 36. Addition of turbulence results in a higher 
temperature in both the H II region, and the partially-ionized zone. In the H II region, 
most of the temperature increase occurs in the neutral He/ionized hydrogen region, and 
is a consequence of the disappearance of that region; that is, the transition from He"*" to 
neutral helium moves closer to the hydrogen ionization front. There is an overall increase in 
the level of ionization in the H II region, including a decrease in the Fe"*"^ column. This is 
probably a consequence of the fact that ions are better able to absorb continuum photons, 
and so the energy of the continuum is released into the gas. In the partially-ionized zone, the 
temperature increase is a consequence of the loss of the synergy described in §4. Once the 
temperature becomes high enough that metal ions can be coUisionally excited, the amount 
of energy in the continuum decreases, the excitation temperature decreases, and the energy 
previously locked up in the excited hydrogen and continua is released to heat the gas. 

Turning to the metallicity, we find that the temperature is lower for higher metallicities 
only for the highest ionization regions of the gas. In this region, the optical depth is relatively 
low and the photons escape, cooling the gas. In the intermediate-ionization regions, and the 
partially-ionized zone, the temperature is higher for higher metallicity. In this case, the 
region where hydrogen is ionized and helium is neutral becomes larger as the metallicity 
increases. Practically all metal-ion columns increase along with metallicity, although with 
different slopes. Photoionization of a metal ion releases more energy than a hydrogen ion, 
and in regions where the gas is starting to have appreciable optical depth, those photons may 
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not be able to readily escape to cool the gas. In the partially-ionized region, the temperature 
increases because the availability of more metal ions to excite means that not so much energy 
goes into the continuum and excited-state hydrogen, so it is released into the gas. 

Next we turn to the lines. The behavior of the lines as a function of turbulence has 



already been investigated by iBottorff et al. I (120001 ): we look at only a few lines that behave 



differently for the PHL 1811 and K97 continuum for this particular choice of ionization 
parameter. As discussed in §4.1.3, Mg II is weak in PHL 1811 because the temperature is 
low. As seen in the top left panel of Fig. 37, Mg II increases more quickly with increasing 
turbulence in gas illuminated by the PHL 1811 continuum compared with the K97 continuum 
as a result of the increasing temperature in the line emitting gas. This is also true for the 
coUisionally excited line Si II A1814. In contrast, the Si II A1263, which is predominately 
excited by continuum pumping, increases at the same rate in both cases. Next, we turn to 
the semiforbidden lines Si III] A1892 and C III] A1909. The transitions producing these lines 
have low gf values, so they do not become optically thick, and should not be influenced by 
lower opacity of the turbulent gas. Indeed, for the gas illuminated by the K97 continuum, 
there is very little dependence on turbulence. However, there is fairly strong increase in the 
flux of these coUisionally-excited lines in gas illuminated by the PHL 1811 continuum that 
is a consequence of the increase in temperature. 

The results when the metallicity is varied are somewhat different (Fig. 37, right). Gen- 
erally speaking, the generally prominent high-ionization metal lines (e.g., C IV) do not 
vary much with metallicity; their emission is primarily a function of gas temperature. The 
intermediate- and low-ionization metal lines increase with metallicity, generally speaking, 
although with a range of slopes. Mg II increases more quickly with metallicity for the 
gas illuminated by the PHL 1811 continuum compared with the gas illuminated by the K97 
continuum; this is probably a consequence of the increase in gas temperature in the partially- 
ionized zone as the metallicity is increases. Both Si III] and C III] increase with metallicity, 
although with different slopes, again pointing to a lack of robustness as a density indicator. 
Also, the increase is more rapid for the K97 continuum for many of the intermediate ioniza- 
tion lines. This seems to be a consequence of changes in the ionization structure of the gas 
as well as the change in temperature. Finally, we examine the "spike" and "gap" Fe II bands 



( iBaldwin et al. II2004I ). We see that both increase more rapidly for the PHL 1811 continuum, 
but with the "spike" (resonance lines and low-lying levels) increasing more rapidly than the 
"gap" (high-excitation Fe II). Recall that the model spike/gap ratio is much larger than 
observed; thus a high metallicity only exacerbates the problem. 



To summarize this section, we have investigated the effects of turbulence and metallicity 
changes on the gas properties and emission. In order to obtain a qualitative understanding. 
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we have limited ourselves to a single representative ionization parameter and density. We 
find that as the metallicity and turbulence increase, the temperature increases and the energy 
in the continuum decreases, and this generally causes an increase in the emission line fluxes. 
Naturally, lines that are collisionally excited are more sensitive to the temperature increase. 
Weak collisionally-excited line emission is a characteristic of the PHL 1811 spectrum, so 
these results may mean that a high degree of turbulence and high metallicity are not present 
in this object. 
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Table 1. Observing log 



Spectrometer 


Date 


Exposure 


Bandpass 


Slit Width 


Resolution 






(seconds) 








HST STIS FUV MAMA (G140L) 


2001-12-03 


2640 


1119-1715 A 


0.2" 


I.2A 


HST STIS NUV MAMA (G230L) 


2001-12-03 


1667 


1576-3159 A 


0.2" 


5.2A 


HST STIS Optical CCD (G430L) 


2001-12-03 


240 


2894-5702 A 


0.2" 


5.5A 


KPNO 2.1m Goldcam 


2001-07-23 


2700 


6080-9 160A 


1.9" 


4A 


KPNO 2.1m Goldcam 


2001-10-24 


900 


3938-7597 A 


1.9" 


5A 


KPNO 2.1m Goldcam 


2001-10-25 


900 


5466-8532 A 


1.9" 


4A 



Table 2. Emission- line properties 



Emission Line 


Line Model"* 


Measured Wavelength"''^ 


Lab. Wavelength •'^ 


Flux 


Equivalent Width 


Velocity Width" 






f I \ 
(A) 


(A) 


/ir>— 14 —2 — 1\ 

(10 erg cm s ) 


(A) 


(kms ) 


He I 


L (UL) 


7067.2 (f^ ) 


7067.2 


< 0.048 


< 0.17 


1752(f) 




L 


6564.12 ± 0.04 


6564.7 


100.6 ± 0.2 


302 


1752 ± 6 




L 


4864.2 ± 0.09 


4862.7 


29.4 ± 0.2 


50 


1943 ± 19 




L 


4347.1 ± 0.3 


4341.7 


9.5 ±0.3 


12 


1716 ± 60 


Na Id 


L(2,thick) 




5893 58 


1 95 + 06 


4.9 


1590 + 125 


Na Id 


L(2, thick, +He I) 


5894.5 ± 1.1 


5893.58 


1.61 ± 0.12 


4.1 


1365 ± 180 


He I 


L(+Na ID; UL) 


5879.2 ± 0.7 


5877.45 


0.63 


1.6 


1752(f) 


Ca IIh 


G 


3936.6 ± 0.4 


3934.8 


1.30 ± 0.12 


1.3 


993 ± 82 


Ca IIk 


G 


3971.4 


3969.6 


1.02 ± 0.10 


1.1 


993 


Mg II 


G(2,thick) 


2800.3 ± 0.4 


2800.0 


28.2 ± 0.5 


12.9 


2550 ± 110 


ClII] 


integrated flux UL 


1915.8 


1908.7 


< 7.3 


< 1.7 




AlIII 


G(2,thick) 


1850.5 ± 0.8 


1858.8 


8.3 ± 0.4 


1.9 


3210 ± 340 


Fe II19I 




1783.18 ± 0.25 


1786.3 


3.27 ± 0.13 


0.74 


1140 ± 120 


Si II 


T ("UL") 


1815.02 (f) 


1815.02 


< 6.1 


< 1.4 




He II 




1640.41 (f) 


1640.41 


< 5.2 


< 1.11 




CiV 


T 


1549.5 


1549.5 


32.0 ± 0.7 


6.6 




unidentified'' 


G 


1395.2 ± 1.7 




28.0 ±3.0 


5.0 


8470 ± 710 


ClI 


T (UL) 


1335.31 (f) 


1335.31 


< 7.9 


< 1.3 




Si II 


T 


1302.0 ± 0.95 


1307.62 


14.9 ±0.8 


2.4 




Sill 


T 


1259.9 ±1.8 


1263.32 


52.0 ±1.1 


7.9 




nV 


T 


1236.3 ± 0.4 


1240.15 


60.5 ±1.3 


9.0 




Lya 


T 


1216.0 ±0.5 


1215.67 


41.5 ±1.2 


6.0 




Si II 


T 


1193.6 ±0.6 


1194.12 


25.6 ±1.0 


3.6 




cm* 


T 


1175.4 ±0.5 


1175.67 


13.8 ±0.9 


1.9 




Si III* 


G(l) 


1110.9 ±0.5 


1111.59 


1.5 ±0.3 


0.20 


1160 ± 230 


Si III* 


G(6) 


1110.5 ±0.3 


1111.59 


1.7 ±0.03 


0.22 


690 ± 130 


Nil 


G(l) 


1084.2 ± 0.2 


1085.12 


0.90 ± 0.20 


0.11 


460 ±110 


Nil 


G(6) 


1084.2 ± 0.2 


1085.12 


0.95 ±0.02 


0.12 


370 ± 150 



''Codes for line models are as follows: L=Lorentzian profile; UL=upper limit (Ax^ = 6.63); "thick" means that we set multiplet ratios equal to one; 
G=Gaussian line profile; T=C IV template was used for the profile; a number gives the number of components in the fit, usually the number of multiplets, 
or "1" if fitted with a single line. 

''For multiplets, we list the gf weighted values, unless marked "thick" in column 2, in which case we give the average wavelength. 



"Vacuum wavelengths. 
"^Uncorrected for detector resolution. 

•^(f) means that the parameter was fixed at the value given. 
'Using the Vcron-Cetty, Joly, & Veron Fe II template subtraction, 
s Using the Boroson & Green Fe II template subtraction. 

''This 1400A feature is Ukely to be a blend of Si IV and O IV]; since the identification is ambiguous, it is labeled as unidentified. 
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Table 3. LOG Model Results 



Emission Line 


Composite 


Spectra Equivalent Widths^ 


K97 LOG'' 


PHL 1811 LOC'= 


PHL 1811 Observed 




Francis 


Zheng 


; Brotherton 


Vanden Berk 








Lya+N V feature 


52 




87 




94.6 






Lya 




85 




92.9 


90.2 


26.5 


6.0 


nV 




10 




1.11 


4.4 


0.06 


9.0 


1400 A feature 


10 


7.3 


8 


8.1 


5.9 


0.3 


5.0 


cIV 


37 


59 


33 


23.8 


47.2 


0.9 


6.6 


He II 


12<: 


3.9 


7.0== 


0.51 


12.7 


0.24 


< 1.1 


1900 A feature 


22 




17 




14.2 


3.0 




Al III 




3.5 




0.40 


1.3 


0.14 


1.9 


Si III] 




3.5 




o.ie"* 


6.4 


1.1 




ClII] 




17.0 




21. 2== 


6.4 


1.7 


< 1.7 


Mg II 


50 


64 


34 


32.3 


107 


15.0 


12.9 



Note. — All values are equivalent widths in A. 

'^Re ferences for composite spectra are iFrancis et al~l lll99ll) . lzhene et aTl lll997h . lBrotherton et al. I 1 I2OO1I) . and lVanden Berk et al.1 
1I2OOJ) . 

''A covering fraction of 0.25 was assumed. 
'^He II A1640+O III] A1663. 
'^Includes also part of Fe III UV 34. 

"Includes also part of Fe III UV 34, UV 68, and UV 61. 
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Fig. 1. — Merged UV and optical spectra. Prominent lines expected in active galaxies are 
labeled, along with a number of metastable lines. Length of tick marks is proportional to gf 
for individual multiplets. Absorption lines, almost all of which are i dentified as originati ng 
in our Galaxy or in intervening absorption systems, are identified in iJenkins et al. I (120031 ). 
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5500 



6000 6500 7000 
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Fig. 2. — Optical spectra near the H/3 (a) and Ha regions (b). The top hnes show 
the original spectra. The offset middle lines sh ow the Fe II subtracted spe c tra, w here 
the blue and red li nes sh ow the results using the IVeron-Cetty. Joly. fc Veron I (120041 ) and 
( 19921) iron eraission templates, respectively (bottom lines). Note that 



Boroson fc Green 



the I Veron- Getty. Joly. fc Veron I (12004J ) template subtraction yields a more symmetric Ha 
line. 
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IRAS 13224-3809 



5800 6000 6200 

Wavelength (Angstroms) 




3800 3900 4000 4100 

Wavelength (Angstroms) 



4200 



Fig. 3. — Wc observe unusual very low-ionization line emission from Na ID (a) and Ca II 
K&H (b) in PHL 1811, but wc observe no convincing evidence for He I. For comparison, the 
offset spectrum from the narrow-line Seyfert 1 galaxy IRAS 13224—3809 is shown. 
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Fig. 4. — A comparison of the Ca II H&K (grey lines) and H/? (black line) profiles in velocity 
space shows that the Ca II lines are distinctly narrower. 
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Fig. 5. — UV Fe II in PHL 1811 and other quasars. The spectra shown were fitted between 
2650 and 3090 A with a model that included two Gaussians to model the Mg II emission 
at 2800A and the Fe II template plus a line to model the continuum. Then the model 
was extrapolated toward shorter wavelengths, (a) The extrapolated model leaves excesses 
between 2300 and 2500 A, originating in high-excitation Fe II, and shortward of 2200 A, 
originating in Fe III. PHL 1811 falls into this category, {h) the extrapolated model fits these 
spectra well. 
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Fig. 6. — Top: The upper level energies for Fe II emission lines, color coded according to 
their Aij value. Middle: the XFx spectrum of PHL 1811. Bottom: Same as the top panel, 
for Fe III. 
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Fig. 7. — The near-UV spectrum of PHL 1811. Positions of permitted, semiforbidden and 
forbidden lines commonly seen in AGN spectra are marked. All of these lines are present 
in the I Zw 1 spectrum, but forbidden and semiforbidden lines are absent in PHL 1811. A 
feature near 1915A is near the expected position of C III] A1909. The inset shows that the 
centroid of this feature is redward of the expected position of C III] , and it is more likely to 
be the 1914. lA component of Fe III UV34 (see text for details). 
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Fig. 8. — A comparison of the PHL 1811 continu um (black line) with the Large Bright 
Quasar Survey composite spectrum (hght grey hne; iFrancis et al. Ill99ll ). Notable features 
are the lack of a break in the PHL 1811 continuum toward short wavelengths in the far UV, 
and the large excess iron emission in the near UV. Insets show a comparison of the C IV 
lines and the H/3 region. 
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Fig. 9. — A comparison of the far UV spectra of AGN. Left: PHL 1811 and similar high- 
luminosity narrow-line quasars RX J0134.2 — 4258, and PHL 1 092. Right: Comparison 
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ward of ~ 1450A, very weak C IV emission, and a severely blended low equivalent width 
Lya/N V feature. 
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Fig. 10. — The region of the spectrum near C IV. The rest positions of the C IV doublet 
hues are marked by dotted hues. Also marked are the positions of plausible weak Galactic 
absorption lines, a nd an absorpt i on lin e from C IV in the z{abs) = 0.17651 intervening 



system reported in I Jenkins et al 



(120031 ). The template model for the C IV profile is shown; 
the dashed lines show th e doublet components, der ived from the sum profile (solid hne) using 
the method described in iLeighly Sz Moore 1j2004h. The dot-dashed hne shows the template 



obtained for the NLSl IRAS 13224-3809 jLeighlv fc Moore Il2004[ l. 
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Fig. 11. — Modeling the Lya/N V feature. Top: The far UV spectrum with the absorption 
hues (grey) modeled out. Middle: The fit using five templates corresponding to Lya, N V, 
Si II A1308, Si II A1263, and Si II A1194, fixed at their rest wavelengths shown by the vertical 
grey lines. Bottom: The fit using the five templates above, and one for C III* A1176 as well, 
in which the wavelengths of the templates were free parameters. The detected metastable 
Si III* A1112 is marked as well. See §3.3.2 of the text for details. 
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Fig. 12. — Comparison of spectral energy distributions. The dashed hnes (dotted hnes 
in the unobservable EUV) show SEDs derived fror n coordinated observations of PHL 1811 
and the narrow- hne Seyfert 1 galaxy RE 1034+39 JCasebeer. Leighlv k Bar"^l2006h . The 
solid black line shows the K97 continuum intended to represent the SED of a typical 



AGN ( iKorista et al. Ill997l). The thick grey lines sh ow the semiempirical SEDs used for 



Cloudy simulations (ICasebeer. Leighly fc Baron II2006I ). labeled by temperature of the big- 
blue-bump rollover, that correspond most closely to the observed SEDs {kTcut = 10 eV 
for PHL 1811 and kT^ut = 230 eV for RE 1034+39, respectively), and to the K97 SED 
{kTcut = 80 eV). The SEDs used for the Cloudy simulations correspond reasonably well with 
the observed SEDs and K97 SED, at least in terms of the position of the peak of emission. 
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Fig. 13. — Gross properties of the SEDs used for simulations. Information for each SED is 
plotted as a function of the SED, where the CLB continua are coded by the b ig blue bump 
cutoff temperature (leftmost 32 points; see ICasebeer. Leighly fc Baron I l2006l . for details), 
and the two points on the r ight-hand side are for the observed PHL 1811 continuum, and 



the AGN continuum used by lKorista et al. I (119971 ). The top panel shows the aox, the middle 
panel shows the continuum flux at 2500A for a photoionizing flux of 10^° photons s~^ cm~^, 
and the lower panel shows the Compton temperature. The PHL 1811 continuum is char- 
acterized by a very steep aox, and requires a relatively bright optical and UV in order to 
produce the same photoionizing flux as the other continua. This will reduce the line equiva- 
lent widths. The soft spectral energy distribution produces a very low Compton temperature 
of around 5 x lO'^ K. 
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Fig. 14. — The electron temperature as a function of depth for our range of spectral energy 
distributions and the fiducial ionization parameter logf/ = —1.5 and density logn = 11. 
The thin lines show the CLB continua with light to dark showing the progression of soft 
SEDS {kT^ut = 10 eV) to hard ones {kT^ut = 320 eV). The sohd thick line and the dash-dot 
line shows the result for the PHL 1811 continuum and the K97 continuum, respectively. The 
electron temperature is lower for the PHL 1811 continuum than any other continuum. 
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Fig. 15. — The hydrogen ionization fraction in the vicinity of the hydrogen ionization front 
and into the partially-ionized zone for the range of spectral energy distributions and the 
fiducial ionization parameter logU = —1.5 and density logn = 11.0. The lines have the 
same meaning as in Fig. 14. For the CLB continua, the fraction of ionized hydrogen at 
a particular depth decreases as the continuum becomes softer {kTcut decreases), but then 
increases again for the softest continua. Likewise, the hydrogen ionization fraction is higher 
for the PHL 1811 continua. The resulting total ionized hydrogen column is 25% greater for 
the PHL 1811 continuum than for the K97 continuum. 
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Fig. 16. — The rate of ionization of hydrogen from n — 2 from coUisional ionization (left 
panel) and photoionization (right panel). The upper group of lines in the photoionization 
plot shows the rate from the direct continuum, and the lower group of lines shows the rate 
from the diffuse emission. The lines and shading have the same meaning as in Fig. 14. The 
rates are dominated by photoionization from the direct continuum, and are very large for 
PHL 1811 and other soft SEDs. 



-79- 



Cl 




Depth Depth 

Fig. 17. — Left: The density of hydrogen in n = 2 as a function of depth for the fiducial 
parameters (logf/ = —1.5 and logn = 11). The hues have the same meaning as in Fig. 14. 
In the partially ionized zone, the density of excited hydrogen is larger in gas illuminated by 
hard SEDs than in gas illuminated by soft SEDs, but recovers for PHL 1811. Right: The 
ratio of the excitation temperature to the electron temperature for the fiducial parameters. 
The ratio is near one for the hard SEDs, implying that the number of hydrogen inn — 2 are 
consistent with that expected from the electron temperature of the gas, and therefore the 
excitation mechanism is collisions. In contrast, the excitation of PHL 1811 is higher than 
predicted by the rather low temperature. 
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Fig. 18. — The ratio of the density of hydrogen in n = 5 from the Cloudy output to the density 
expected based on the density in n = 1 and the electron temperature for the fiducial gas 
parameters (logf/ = —1.5 and logn = 11); this parameter is proportional to the departure 
coefficient for n = 5. The lines have the same meaning as in Fig. 14. The departure coefficient 
is very large for PHL 1811. 
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Fig. 19. — Integrated diffuse continuum times the depth (in ergs^^cm"^) as a function of 
depth for the fiducial parameters (logC/ = —1.5 and logn = 11). The hues have the same 
meaning as in Fig. 14. The energy in the diffuse continuum is greater in gas illuminated by 
the hard continuum compared with the gas illuminated by the softer continua, but is even 
higher in PHL 1811. The energy in the continuum is large in objects with hard SEDs because 
of the greater amount of energy injected into that region by the soft X-ray photoionization. 
The energy in the continuum is large in PHL 1811 because the gas is so cool that usual 
channels of cooling (e.g., excitation of metal ions) are less important, and the gas has no 
options besides emission of continuum. The gas illuminated by the hard continuum resembles 
gas with a high density in that cooling options are limited. 
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Fig. 20. — A comparison of the predicted ion and emission- line properties for the PHL 1811 
continuum and the K97 continuum for the fiducial gas parameters (log U — —1.5 and logn = 
11.0). In both cases, the column density is 10^^-^, which is sufficiently large to be very 
optically thick to the hydrogen continuum. The embedded small panels show the same 
information as the larger panels on an expanded scale. The top plot shows the ratio of 
the predicted ion columns for PHL 1811 to K97 as a function of the ionization potential. 
For the lowest ionization potentials, the columns are nearly the same. For the highest 
ionization potentials, the column is much lower for PHL 1811 continuum, as it is too soft 
to create or excite highly-ionized ions. Intermediate-ionization ions fill the region above the 
hydrogen ionization front, and thus their levels are elevated in PHL 1811. The bottom panel 
shows the ratio of resulting common emission lines. The disperson for a lines emitted by 
a particular ion arises from the fact that lines that are coUisionally excited are generally 
weaker for PHL 1811 because the temperature is lower, whereas lines that are excited by 
continuum fluorescence are brighter because the UV continuum is relatively brighter for the 
same ionization parameter and density in PHL 1811 compared with K97. 
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Fig. 21. — A comparison of the line emissivity and ratios for the PHL 1811 continuum and 
the K97 continuum for Si II A1814 (grey hues) and Si II A1194 (black lines). The top panel 
shows the line emissivity multiplied by the depth as a function of the depth. For the 1814A 
line, the emissivity of PHL 1811 is smaller than that of K97 particularly in the partially- 
ionized zone where Si"*" dominates (depths >~ 10^^ cm). In contrast, the emissivity of the 
II94A line is much higher for PHL 1811. The lower panel shows the ratios from the gas 
illuminated by the PHL 1811 continuum to that of the gas illuminated by the K97 continuum. 
The solid lines show the ratios computed from the Cloudy predictions of the line emissivity 
as a function of depth, while the dashed lines show the ratios computed using the electron 
temperatures and the Boltzmann excitation equation. The two ratios match well for the 
collisionally-excited line at 1814A; deviations at low depths are a consequence of different 
fractional abundances of Si"*" ions. But the observed ratio is much larger than the Boltzmann 
equation ratio for the 11 94 A hne. This is because the 1194A line is predominately excited 
by continuum fluorescence, and much brighter in PHL 1811 due to its relatively brighter UV 
continuum for the same ionization parameter and density. 
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Fig. 22. — A comparison of the line equivalent widths of C IV for the PHL 1811 continuum 
(left) and the K97 continuum (right) as a function of ionization parameter and density. C IV 
is predicted to be weak in PHL 1811 because of the low column of C^^ ions and the lower 
temperature in the emitting gas of this coUisionally excited line. 
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Fig. 23. — A comparison of the line equivalent widths of Si III] for the PHL 1811 continuum 
(left) and the K97 continuum (right) as a function of ionization parameter and density. This 
semiforbidden line has a critical density of 3 x 10^^ cm~^, so a decline toward higher densities 
in the K97 simulation is expected. In contrast, in the PHL 1811 continuum simulations, the 
onset of the decline toward higher densities occurs at much lower densities, around 10^° cm^"'. 
This implies that for gas illuminated by soft continua, the usual density indicators are not 
reliable. 
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Fig. 24. — Properties of Fe II using the Cloudy 371-level atom for the fiducial parameter 
(logf/ = —1.5 and logra = 11). Top: The ratio of the Fe II emission in broad bands for the 
range of spectral energy distributions compared to that of PHL 1811. Solid circles: 1000- 
2000A, solid stars: 2000-3000A, x's: 4000-6000A, open circles: 7800A-3 microns. Second 
from top: Ratio of Fe II flux in the 2200-3050A band to the Mg II flux. The large ratio 
for PHL 1811 is a combination of higher Fe II flux, a consequence of stronger continuum 
pumping, and lower Mg II flux, a consequence of lower electron temperature. Third from 
the top: The ratio of the Fe II emission in the UV (2000-3000A) to the Fe II emission in 
the optical (4000-6000A). This ratio is high in PHL 1811 because the continuum pumping 
enhances the UV emission, and because of Fe II emission in the H+-neutral He region. 



Bottom: The "spike/gap" ratio, defined by iBaldwin et al. I (l2004l ). This ratio is essentially 
the emission from the low-lying, collisionally-excited levels that produce the spikes, to the 
emission from high levels. The high-level emission is enhanced in PHL 1811 due to continuum 
pumping. 
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Fig. 25. — Properties of level number 295 (10.85 eV) as a function of depth for a range 
of spectral energy distributions for the fiducial parameters (logf/ = —1.5, logn = 11, and 
log Nh = 24.5). The lines have the same meaning as in Fig. 14. Left: The departure co- 
efficient shows that level 295 is much more highly populated (by a factor of almost 100 in 
the partially-ionized zone) in the gas illuminated by the PHL 1811 spectral energy distribu- 
tion than expected from the local electron temperature. Right: The density of Fe+ ions in 
level 295 shows that the departure coefficient is high for gas in PHL 1811 both because the 
temperature is low (Fig. 14) and because the number of atoms with electrons in this level is 
high. 
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Fig. 26. — The emissivity times depth as a function of depth for two regions of the Fe II 
spectrum for the fiducial parameters (logf/ = —1.5, logn = 11, and logNn = 24.5). The 
hues have the same meaning as in Fig. 14. Note that the emissivity has the opacity factored 
in; the value is the amount reaching the surface. The left side shows the "spike" region of 
the spe ctrum (2312-2328 an d 2565-2665A) and the right side shows the "gap" region (2462- 
2530A; [Baldwin et aini2004h . In both cases, the "hook" between ~ 3 x 10^° and ~ 10" cm 
for PHL 1811 originates in the H"'"-neutral He region. The spike is predominately resonance 
transitions, and the emission in the partially-ionized zone is approximately proportional to 
the temperature of the gas (although somewhat elevated for the softest continua). The gap 
is predominately emission from higher levels; these levels are pumped by the continuum and 
Lja when the SED is soft, and they are higher in the partially-ionized zone. 
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Fig. 27.— The log of the ratio of the Fe II emission in 2200-3050A band to the Mg II 
emission as a function of ionization parameter and density. In each case, the column density 
is proportional to the ionization parameter such that log Nh + log C/ = 26. The ratio is 
consistently higher for PHL 1811 (left) compared with the K97 continuum (right). 
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Fig. 28.— The "spike/gap" ratio for the PHL 1811 (left) and K97 (right) SEDs as a function 
of ionization parameter and density. The ratio for the PHL 1811 continuum is consistently 
lower than for the K 97 continuum, although still much larger than the observed value (0.7; 
Baldwin et al. II2004I ). As discussed in the text, the lower ratio originates from higher emis- 
sion from high levels that comprise the "gap" emission in gas illuminated by the PHL 1811 
continuum. 
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Fig. 29. — Contours for Cloudy simulations of Ha, the upper limit for He I A5875, and their 
ratio. Only regions in which the log of the ratio of the simulated to observed is between —1 
and 1 (i.e., 1 dex from a perfect match) are shaded. Four parameters {kTcut, D, \og{density) 
and \og{U) define the parameter space of the simulations, where D is ratio of the depth to 
the depth of the hydrogen ionization front. Shown are contours in planes defined by two of 
the four parameters: kTcut and D on the left side, and \og{density) and log(f/) on the right 
side. These planes intersect at the point kTcut = 10, log{density) = 11.5, log(6'^) = —2.0, 
and D = 10. 
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Fig. 30. — Results from Cloudy simulations in which the column density is chosen to max- 
imize the ratio of the He'*' and H^" columns. Left: The log of the ratio of the transmitted 
continuum just below the He I edge (23.9 eV) to the transmitted continuum just below the 
N II edge (28.8 eV). A plateau exists upon which the He I edge is very prominent in the 
transmitted spectrum. Right: A representative example continuum for log(L'^) = —1.5 and 
log(n) = 11.5. The dashed black line shows the PHL 1811 continuum. The light grey line 
shows the continuum after it is absorbed by the gas. The dark grey line shows the absorbed 
plus diffuse continuum, or the transmitted continuum. The inset shows a magnified view 
of the region around the Lyman and He I edges. The diffuse continuum fills in the Lyman 
edge. 
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Fig. 31. — The results of filtering the continuum on the He I A5786 emission for the case 
of the PHL 1811 continuum. The left column, labeled (a.), shows the He I emission with 
respect to the upper limit when the emitting gas is illuminated directly by the continuum 
(no filtering). A covering fraction of 10% is assumed. Except for low column densities, the 
predicted He I exceeds the upper limit by a large factor. The right column, labeled (b.), 
shows the He I emission with respect to the upper limit when the emitted gas is illuminated 
by the transmitted continuum shown in Fig. 30 (the filtered continuum). The predicted He I 
emission is lower than the upper limit for low to moderate ionization parameters. In both 
columns, the top panel shows the case where log(n) = 10.5, and the bottom panel shows the 
case where \og{U) = —2.75. Note that the ionization parameter shown is the one appropriate 
for the unfiltered continuum in both cases. 



-94- 




0.01 0.10 1.00 10.00 0.01 0.10 1.00 10.00 
Depth Depth 




-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 
log(simulated/observed) 



Fig. 32. — The corresponding results of filtering the continuum on the Ha emission for 
the same case as shown in Fig. 31. In both columns, the top panel shows the case where 
log(n) = 10.5, and the bottom panel shows the case where \og{U) ~ —2.75. Again, a covering 
fraction of 10% is assumed. Note that the ionization parameter shown is the one appropriate 
for the unfiltered continuum in both cases. There is httle difference between the non-filtered 
case (left) and the filtered case (right). This is expected because the Lyman continuum is not 
much altered by the transmission, as shown in Fig. 30. Ha is in the vicinity of the observed 
value (although lower by a factor of two) at high column densities and moderate densities. 
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Fig. 33. — The spectral energy distributions used in Appendix A.l. The SEDs are based on 
the PHL 1811 and are described by two parameters: aox-, and the break between the UV and 
the X-ray. All SEDs are shown for the same value of the photoionizing flux. These SEDS 
span parameter space described by the strength of the extreme UV and the strength of the 
X-ray between the PHL 1811 (thick sohd line) and the K97 (thick dash-dot hne). 
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Fig. 34. — Contours of the log of the hne flux as a function of the parameters describing 
the SEDs shown in Fig. 33, aox and the EUV cutoff, evaluated for the fiducial parameters 
(logf/ = —1.5, logn = 11.0, and logA'^^ = 24.5). A variety of lines are shown, including 
lines emitted in the H II region (top row), the partially-ionized zone (middle row), and UV 
Fe II lines (bottom row). In each graph the shaded contour interval is 0.05. The behavior 
of high-ionization lines is governed by the ionization potential of the lines. CoUisionally 
excited intermediate- and low-ionization lines such as C III], Si III], Mg II, and Si II A1814 
depend primarily on the tempera t ure. P umped lines including Si II A1194, and Fe II lines 
UV 191, and the ISaldwin et aTl J2004h "spike" and "gap" lines depend on the UV flux, 
which is stronger for a given photoionizing flux for lower values of the EUV cutoff, and the 
temperature. 
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Fig. 35. — Properties of gas illuminated by SEDs shown in Fig. 33, evaluated at a depth 
of ten times the depth of the hydrogen ionization front. The electron temperature (left), 
the hydrogen ionization fraction (middle) and the ratio of the excitation temperature to the 
electron temperature for n = 2 are shown as functions of the parameters describing the 
SEDs, and the EUV cutoff. 
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Fig. 36. — The electron temperature as a function of the ratio of the depth to the depth of 
the hydrogen ionization front for the fiducial gas parameters (logC/ = —1.5, logn = 11, and 
log Nh — 24.5) as microturbulcnce (left) and metallicity (right) are varied. The solid line 
shows the result for the PHL 1811 continuum with Vturb = and Z/Zq = 1. Light to dark 
shows the results as the parameters are increased. Increase of either turbulence or metallicity 
results in an increase in temperature in the partially-ionized zone to values comparable to 
those obtained for typical AGN continua (Fig. 14). 
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Fig. 37. — The ratio of lines as a function of turbulence (left) and metallicity (right) to the 
values obtained for models in which Vturh = and Z/Zq = 1, respectively, evaluated for 
the fiducial parameters {\ogU = —1.5, logn = 11.0, and logNn = 24.5). The solid and he 
dashed lines show the case for the PHL 1811 and the K97 SEDs. More rapid increases in 
the fine ratios for the PHL 1811 SED compared with the K97 SED are a consequence of 
temperature increases with turbulence and metallicity. 



